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Abstract. Investigating the stability of nonlinear waves often leads to lin- 
ear or nonlinear eigenvalue problems for differential operators on unbounded 
domains. 

In this paper we propose to detect and approximate the point spectra of 
such operators (and the associated eigenfunctions) via contour integrals of so- 
lutions to resolvent equations. The approach is based on Keldysh' theorem 
and extends a recent method for matrices depending analytically on the eigen- 
value parameter. We show that errors are well- controlled under very general 
assumptions when the resolvent equations are solved via boundary value prob- 
lems on finite domains. Two applications are presented: an analytical study 
of Schrodinger operators on the real line as well as on bounded intervals and 
a numerical study of the FitzHugh-Nagumo system. 

We also relate the contour method to the well-known Evans function and 
show that our approach provides an alternative to evaluating and computing 
its zeroes. 
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1. Introduction 



Studies of the analytic Evans function ([AGJ, PW]) have found numerous apph- 
cations to the stabihty problem for nonlinear waves. The function has proved to be 
an invaluable tool for locating point spectra of differential operators that are defined 
on the real line or, more generally, on cylindrical domains. In particular, analyzing 
the behavior of the Evans function near zero or near infinity provides one major 
step in studying point spectra rigorously. We refer to [Sa], [KKS], [GLM],[OS], [DL] 
for a variety of applications and extensions of the concept to other settings. 

There have been quite a few approaches to computing the Evans function numer- 
ically and then determine from its zeros the unknown eigenvalues, see [Br], [BDG], 
[BZ], [HZ], [HSZ], [MN], [LMT], [SE]. The standard definition of the Evans func- 
tion (see [Sa]) involves a determinant of vectors which depend analytically on the 
eigenvalue parameter and which are determined as initial values of exponentially 
decaying solutions on both semi-axes. This leads to the problem of integrating 
stiff ODEs while keeping the analyticity with respect to the parameter. For higher 
dimensions the evaluation of a determinant may also lead to instabilities unless a 
proper scaling is employed, see for instance [GLZ, Thm.4.15] for a discussion of the 
appropriate scaling of the Evans function. The techniques proposed to solve these 
problems utilize exterior products [Br],[BZ],[AB], [BDG] or solve Kato's matrix 
differential equation [K, Sec. II. 4. 2], see [HSZ]. 

In this paper we pursue an alternative road that avoids the intermediate stage of 
computing the Evans function. Rather we propose to solve the original analytical 
eigenvalue problem in the whole space via well-posed boundary value problems on 
finite domains while preserving analyticity. For the latter problem we propose a 
numerical method using contour integrals of solutions to resolvent equations. The 
approach is based on Keldysh' theorem and extends a recent method [B] for non- 
linear eigenvalue problems with matrices. Om- goal is to determine all eigenvalues 
inside a given contour F c C and to guarantee that, at each stage of approximation, 
approximate equations are as well-conditioned as the original eigenvalue problem. 
The method applies to general nonlinear eigenvalue problems 



where F{\) are Fredholm operators of index that depend analytically on A G fi. 
For certain right-hand sides v and functionals w we require to evaluate the integrals 



where F lies in the resolvent set. The idea first appeared for the matrix case in 
the papers [AK08],[AK09] where it was used in connection with the Smith normal 
form. In [B] we arrived independently at similar expressions (without applying the 
functional w) by using the theorem of Keldysh, see [MM, Ch.l]. It turns out that the 
Keldysh approach allows to locate eigenvalues precisely, to handle multiplicities and 
to simultaneously compute (generalized) eigenvectors, see Section 2.3. Moreover, 
Keldysh' Theorem works for abstract Fredholm operators [MM] and thus allows 
to generalize the whole approach as we will show in Section 2. Note that contour 
integrals are normally used for computing the winding number and thus the number 
of zeros inside the contour (see [Br], [BDG], [BZ] in case of the Evans function), 
whereas our method allows to locate all zeros inside the contour and, in addition, 
to obtain approximate eigenfunctions (Theorems 2.4, 2.8 and equation (2.23)). 



F{X)v = 0, A e O 



(1.1) 




and 




(1.2) 



EIGENVALUES AND BOUNDARY VALUE PROBLEMS 



3 



In Section 3 wc discuss suitable normalizations of the Evans function and relate 
our contour method to such a normalized Evans function (Theorems 3.13, 3.18). 
Then we apply the method to first order differential systems with A-dependent 
matrices. We continue this in Section 4 and show that isolated eigenvalues are well 
approximated when the resolvent equations arc solved on bounded intervals with 
suitable boundary conditions (Theorems 4.7, 4.11). Here we follow [BR] where it 
is shown that the functional analysis of discretization methods [V76] , [V80] applies 
to this situation. In Section 5 we apply the theory to Schrodinger operators on 
the real line and express the various quantities of our approach in terms of Jost 
solutions and Levinson's Theorem (Theorems 5.3, 5.4). Section 6 concludes with 
several numerical experiments for the FitzHugh-Nagumo system which demonstrate 
the robustness of the contour method. In particular, we show how computational 
errors depend on the length of the finite interval, the number of quadrature points 
used for (1.2), and on the rank test used to determine the number of eigenvalues 
inside the contour. 

2. Abstract results for the contour method 

2.1. The abstract setting. We consider nonlinear eigenvalue problems as in (1.1) 

where F : — > £('H, K.) is a holomorphic function on a domain C C with values 
in the space £(7^, fC) of bounded linear operators from some complex Banach space 
into another Banach space JC. We will assume that the operators F{X) are 
Fredholm of index for all A G (this will be written as G H(f}, J"(-H, /C))). As 
usual, cf. [MM, Ch.I], we define the resolvent set and the spectrum of the operator 
pencil F by 

p{F) = {A e n : F(A) invertible}, a{F) = O \ p{F). 

Throughout, we impose the following assumptions. 

Hypothesis 2.1. Assume that F e H(n, J"('H,/C)) and p{F) ^ 0. 

Under these assumptions cr{F) is a discrete subset of fl and the operator valued 
function F~^{-) is meromorphic, see [MM, Thm. 1.3.1]. In the terminology of [GGK, 
Ch.IX] the function F~^{-) is finitely meromorphic. In general, we follow the setting 
in [MM, Ch.I]. In particular, we use dual spaces Ti', K' and denote the dual pairing 
by elements w G'H', u G )C' in two equivalent ways 

w'^'v = {Wjv), vgH, vJv = {u,v), V G K.. (2.1) 

Thus, the dual space; is the space of linear (versus complex conjugate linear) func- 
tionals. As noted in [MM, Sec. 1.1] it is important to avoid the complex adjoint space 
and complex conjugate functionals (even if H and /C are Hilbert spaces), since this 
will turn holomorphic functions into antiholomorphic ones. We prefer the notion 
over w* and w' because it is in accordance with matrix-vector notation. Cor- 
respondingly, we denote the dual of an operator A e £{'H,IC) by e C{IC' ,%'). 
It is defined by 

{A'^w,v) = {A'^w^v = w'''{Av) = {w,Av), w € IC',v gH. 
Any two elements v and w G )C' define an operator e £(/C, H) as follows 
{vw''')u := V {w'^'u) = V {w,u), u&K. (2.2) 
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Because of this definition we omit brackets and simply write vw u. Note that the 
operator vw'^ is of rank one it v 0,w ^ 0. We denote by Af{-) the null-space and 
by TZ{-) the range of a linear operator. 

Let r be a smooth contour in fl surrounding a bounded subdomain int(r) = 
fio C O. Since f2o U F is compact and eigenvalues are isolated [MM, Thm. 1.3.1] 
there are at most finitely many eigenvalues Ai, A2, . . . , A;,^ e ^0, and, in addition. 



they are of finite multiplicity. Our goal is to determine these eigenvalues and good 
approximations of the eigenvectors by computing contour integrals of the type (1.2) 
In the following choose m linearly independent functionals Wj £ H' , j = 1, . . . ,m, 
and £ linearly independent vectors Vk & JC, k = 1, . . . ,£. Let us assume that F{X} 
is invortiblo for all A € F. Then we can solve the following equations for unknown 
vectors yk{^) S H: 

F{X)yk{X)=Vk, A G F, fc = (2.3) 
We introduce an (m x £)-matrix valued function, E{-), on F as follows: 

(tZ;,-,t/fe(A))) Aer. (2.4) 

In addition, we introduce the following (m x i) matrices: 

Do = ^ [ E{\) dX, £)i = -L /■ XE{X) dX. (2.5) 



Remark 2.2. We note that equation (2.3) is the only information used to obtain 
E{X) and thus the matrices £>o and Di. So, replacing (2.3) by an approximate 
equation with good stability properties, the matrices Dq,Di can be computed us- 
ing approximation arguments. In our applications, equation (2.3) is a differential 
equation on the line, while the respective approximation is a differential equation 
on a finite segment with appropriate boundary conditions. Wc will discuss the er- 
rors of this approximation in Sections 2.4 and 4 below. In addition, there are errors 
caused by approximating the contour integrals (2.5) by a quadrature rule, e.g. the 
trapezoid sum. This error depends on the number of quadrature points and on the 
distance of eigenvalues to the contour. For a detailed analysis in case of analytic 
contours F we refer to [B]. In any rate, we may assume in what follows that the 
matrices £^(A), Dq, Di are known to us. O 

Remark 2.3. At first sight it seems unncccesarily general to have different dimen- 
sions m and £ that may also differ from m:. However, it is important for practical 
computations. First of all, >f is generally unknown and thus a quantity to be 
determined. Our approach will only require m,£ ^ x which can be achieved by 
increasing m and £. Second, the number £ determines the number of equations (2.3) 
to be solved and hence the numerical effort. Therefore we want to keep it as small 
as possible. Finally, m can be very large since it gives the number of functionals 
that we can evaluate on the solutions of (2.3). In principle we can assume to know 
yj{X) exactly and hence all values w^t/j(A),ti; e /C'. In this case we may think of 
E{X) being a matrix with infinitely long columns, see Example 4.6. O 

2.2. Simple eigenvalues inside the contour. In this subsection, in addition to 
Hypothesis 2.1, we assume that F has only simple eigenvalues inside the contour 
F. That is, for some x > 0, 

a(i^)nOo = {Ai,...A^}, (2.6) 
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dim J\f{F{\n)) = 1, and there are eigenvectors Vn & of the operator F(A„), the 
eigenvectors Wn € /C' of the dual operator F(A„)' such that, cf. [MM, Def.1.7.1], 

F{\k)vn = 0, wjF(A„) = 0, «;jF'(A„K^O, n=l,...,x. (2.7) 

In the following it will be convenient to normalize Vn,Wn such that 

u;jF'(A„)t;„ = 1, n = l,...,x. (2.8) 

We recall the well-known Keldysh formula, see [MM, Thm.1.6.5]: 

^ W = E tA-^"^" + ^W' A e fio \ {Ai, . . . , A^}, (2.9) 

where the normalization (2.8) is assumed and H{-) is a holomorphic function on a 
neighborhood U of flo with values in £(/C,'H). We also use the notation from (2.2). 

Using (2.9) in (2.4) yields for j = 1, . . . ,m, k = !,...,£ the following formula 
for the entries of the matrix E{X): 

^ 1 

EjkW = Y] T ^{Wj,Vn){Wn,Vk) + {Wj,H{X)Vk{X)). (2.10) 

n=l ^ ~ ^" 

We write this in matrix form by introducing the rectangular matrices Gi G C"*'"^ 
and Gr e C^'^, 

, Gr=[{Wn,Vk)] , (2.11) 

/J — l,n— 1 \ / A;— l,n— 1 

the diagonal matrix A = diag{Ai, . . . , A;^}, and the matrix function 

Then (2.10) has the matrix form 

E{X)=Gi{XI^-A)-^Gj + Ho{X), Xeil, (2.12) 

where Hq € H(W, C"'^). Using Cauchy's Theorem, we evaluate the matrix Dq from 
(2.5) as follows: 

Do = ^.J Gi{XI^ - A)-^GjdX = GiGj. (2.13) 
A similar calculation yields 

£>! = GiAGj. (2.14) 

Our standing assumption in the following will be that both matrices Gi and Gr 
have rank >c. As we noted in Remark 2.3 this requires to have £ ^ x and m ^ >c. 
Equation (2.13) then implies the following formula for the number of eigenvalues 
Afe enclosed by F: 

>f = rankDo. (2.15) 

Thus, cf. Remark 2.2, one can compute the number >c from a rank test of Dq and 
this should be sufficiently robust to approximation arguments. 

Next, wc show how Di can be used to evaluate the actual location of the eigen- 
values An enclosed by F. And this procedure should be robust to approximation 
arguments as well. 



6 



W.-J. BEYN, Y. LATUSHKIN AND J. ROTTMANN-MATTHES 



Let (Ti, . . . , <J^ denote the nonzero singular values of the matrix Do, and introduce 
the diagonal {k x >i) matrix So = diag{cri, .... cTj^}. We use the short form of the 
singular value decomposition of (e.g. [A05, §3.2]) 

= FoSoWo*, e C'"-'^, ^0*1^0 = Wo e C^'", Wo*Wo = 1^. (2.16) 

Note that here one uses adjoint matrices Wq = Wo""" and = Vo^ ■ 
Due to equation (2.13), we have 

Do = Vo^oW; = GiGj. (2.17) 

Since the columns of the matrices G; and Vq span the same subspace, there is a 
nonsingular (x x x) matrix S such that VqS = Gi. Prom this we obtain 

EoWo* = 5G7 and GJ = S-'^'SqW^. 

Using this in equation (2.14) we find 

Di = VoSAS-'^lloW^. 

Multiplying the last equation by Vq from the left and by WqSq ^ from the right, 
yields 

SAS^^ = VoDiWaT^o^ =: D, (2.18) 
which, in turn, implies the desired formula for the Afc's: 

{Ai,...,A^}-a(i?). (2.19) 

We stress again that the spectrum of the matrix D in (2.18), by (2.4), (2.5), can be 
evaluated using only the data from (2.3), and thus can be obtained using approxi- 
mation arguments. We summarize our result in the following theorem. 

Theorem 2.4. Suppose that the operator pencil F G M(fl, J^{'H, K,)) has only simple 
eigenvalues Ai , . . . , A,f inside a simple closed contour T inO. and no eigenvalues on 
the contour. Given a set of linearly independent junctionals wj G H' , j = 1, . . . ,m, 
and vectors Wfe € /C, /c = 1, . . . , with m,£ ^ k, let the matrices Dq, Di € (^mx^ 
be determined by (2.3), (2.4), (2.5) and assume that the matrices Gi,Gr from (2.11) 
have maximum rank. Then the x x x-matrix D = Vq DiWqT^q^ , where Vo,So, Wq 
are given by the singular value decomposition (2.16) of Dq, has ..,X^ as simple 
eigenvalues. 

As soon as A is determined, one can view (2.18), that is, the equation 

DS = Sk, (2.20) 

as an equation for the corresponding to A„ eigenvectors of the matrix £>, which 
are the columns of the matrix S. With S known we can compute Gi = VqS, i.e. 
the values {wj,Vn) are available. From this we determine good approximations of 
eigenvectors n = 1, . . . , x, as follows: 

Step 1: Select vectors Uk & H,k = 1, . . . ,m, that are biorthogonal to the tSj's, 

{wj,Uk) = 5jk, j,k = l,...,m. (2-21) 
Step 2: Determine coefficients /3fc,n, k = 1, . . . ,m, which minimize 

m m 
= J2\i^j^^n-Y.l^k,nUk)\''. (2.22) 
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Due to (2.21) the minimum is attained at l3k,n = {wk,Vn) = {Gi)k,n- 
Step 3: Determine approximate eigenvectors of the operator F(A„) from 

m 

v^^"^'"^ = Y,{Gi)k,nUk, n (2.23) 
fe=i 

If is a Hilbcrt space then we can identify uij = Uj,j = l,...,m, so that 
(2.21) requires these vectors to form an orthonormal system. Introducing the 
subspace = spanjui, . . . ,Sm} we find that ^{13) in (2.22) agrees with \\vn — 
Sa—i /^fe.nWfep up to a constant. Hence ii^pprox jg ^-jj^ host approximation of w„ in 
the subspace X^- If the functions uj are not orthonormal, one can find the best 
approximation of Vn in X^ by first solving for I3j^n the linear system 

rn 

^{uk,Uj)nl3j,n = {wk,Vn)n, k=l,...,m,n=l,. ..,>(, 
and then setting 

m 

^approx^^^.^^^.^ n = l,...,x. (2.24) 

i=i 

2.3. Multiple eigenvalues inside the contour. In this section we show that 
multiple eigenvalues do not produce serious problems with the contour method. 

In Theorem 2.8 we prove that the matrix D from (2.18) inherits the miihiphcity 
structure of the original nonlinear problem. This is analogous to the behavior of 
the Evans function (sec [AGJ]), see Section 3.1 for more details. Let us first recall 
the definition of multiplicity of eigenvalues for nonlinear pencils and the associated 
notion of chains of eigenvectors (see [MM, Sec. 1. 1.6], [V76]). 

Definition 2.5. Let F G H(n, J"(H, /C)) and Aq G a{F). 

(i) A tuple {vq, . . . , Vn-i) € 'H", n > 1, is called a chain of generalized eigenvec- 
tors (CGE) of F at Ao if the polynomial v{\) = J2^Zq{^ — ^oYvj satisfies 

{Fv)^^\Xo)=0, i = 0,...,n-l. 

The order of the chain is the index ro n) satisfying 

(i^t;)W(Ao)=0, i = 0,...,ro-l, (Fu)('-'')(Ao) ^ 0. 

The rank r{vo) of a vector vq € J\f{F{X)),vo 7^ 0, is the maximum order of CGEs 
starting at Vq- 

(a) A canonical system of generalized eigenvectors (CSGE) of F at Aq is a 
system of vectors 

Wj,p 'EH, j = 0,...,fip-l, p=l,...,q, g > 1, 

with the following properties: 

(1) vo,i,. ■ ■ , vo,q form a basis ofAf{F{Xo)), 

(2) the tuple {vo,p, ■ ■ ■ , w^p_i,p) is a CGE of F at Aq for p = 1,. . .,q, 

(3) for p = 1, . . . , g the indices Hp satisfy 

jjLp = max{r(iio) : fo € N{F{\q)) \ span{tio,i. : 1 ^ < p}}. 
(Hi) The numbers jip, p = l,...,q, are called the partial multiplicities, where 
fj-i + ■ ■ ■ + fiq is the algebraic multiplicity and q = dimA/"(-F(Ao)) is the geometric 
multiplicity. The subspace span{v j^p : j = 0, . . . , iJp — l,p = 1, . . . ,q} is called the 
root subspace of the eigenvalue Aq. 
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Such a CSGE always exists and, as with usual Jordan chains for matrices, 
condition (ii)(3) guarantees that chains of highest order are taken first, so that 
Ml ^ At2 ^ • • • ^ holds. 

Next we state the general formula of Keldysh ([MM, Thm.1.6.5]) for a finite 
number of eigenvalues inside a given contour (cf. [B, Cor. 2. 8] for this generaliza- 
tion). 

Theorem 2.6. Let F € EI(f2, J^{H, K)) and let T C p{F) he a simple closed contour 
with hounded interior Q.q in O. Let Oq ncr(F) = {Ai, . . . , A>c} and consider for each 
A„, n = 1, . . . , X, a CSGE denoted hy 

vIp&U, j = 0, p=l,...,qn, n = !,...,>€. 

Then there exist corresponding CSGEs 

wlpGlC', j = 0, ...,/i„,p-l, p=l,...,qn, n = l,...,x, 

for € H(r^, J"(/C', "H')), an open set U with fio U F c C O and a function 
H G m{U,J^{K.,n)) such that 

= EEE(^-^")"'''E'<XZ.-i-. + ^^W' A eZY\a(F). (2.25) 

n—l p—1 j—1 iy—0 

Remark 2.7. The dual CSGEs are uniquely determined by an orthogonality condi- 
tion which generalizes (2.7), see [MM, Thm.1.6.5]. We omit these conditions here 
since they will not be used it in the sequel. O 

With Keldysh' formula (2.25) and Cauchy's formula we repeat the calculations 
that lead to (2.13), (2.14). The result is 

Do = GiGj, Di=GiAGj, (2.26) 

where the matrices Gi,Gr are of size m x >co and £ x hq, respectively with 

(2.27) 



= EE'^"'P' 

n=l p=l 



More explicitly, we have 

{Gr)un,p,.) = {w^„,,-.-i,vj), I ^ o;;;;'/^^_^ - 1, p = i, . . . , g„, n = i, . . . , 

(2.29) 

where one may think of the triples (n, p, v) being ordered lexicographically. More- 
over, it turns out that A is of Jordan normal form 

M \ (Jn,l \ 

A = • . , J„ = 



V 



(K 1 



V 



\ 



K 1 



(2.30) 



(2.31) 
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This immediately leads to the following generalization of Theorem 2.4. 

Theorem 2.8. Let F £ M.{i}, J^{'H, /C)) and let T c p{F) be a simple closed contour 
with bounded interior flo in O. Let Qor\a{F) = {Ai, . . . , Xj^c} and denote the CSGEs 
associated with A„ by 

Vj,p, j = 0, . . .,/i„,p - 1, p=l,...,qn, n=l,...,>c. 

Consider linearly independent elements Wj G T-L' , j = 1, . . . ,m, and Vk € H, k = 
!,...,£, such that m,(. ^ xq, see (2.27). Determine the matrices Dq,Di e C™^^ 
by (2. 3), (2. 4), (2. 5) and assume that the matrices Gi,Gr from (2.28), (2.29) have 
maximum rank. Then the kq x KQ-matrix D = Vq DxWqT,q , where Vb, So) Wq are 
given by the singular value decomposition (2.16) of Do, has Jordan normal form 
(2.30), (2.31) which coincides with the multiplicity structure of the spectrum of the 
nonlinear operator inside T. 

Remarks 2.9. (a) It is well known that the Jordan normal form is not robust to per- 
turbations and, therefore, not a suitable object for numerical computations. There- 
fore, it seems questionable, whether one should compute approximate CSGEs by 
applying formula (2.23) to the columns of Gj in (2.28). Nevertheless, Theorem 2.8 
has some significance. In our case perturbations of the matrices Dq, Di,Vo, Wq, Sq 
may be caused by approximate solutions of the operator equations (2.3), by quad- 
rature errors for the contour integrals (2.5), or by errors of the singular value 
decomposition (2.16)( see Section 2.4 for more details). This leads to a perturbed 
matrix D with well known spectral properties, e.g. the invariant subspaces of D 
belonging to the eigenvalues A„,n = 1, . . . , >c, keep their dimension and perturb 
with the same order, see [StS], and the error of a single eigenvalue A„ grows at 
most with the power l//ii,n of the perturbation, where is the maximum rank 
of eigenvectors belonging to A„ (cf. [K]). 

(b) There are several reasons that may cause a rank defect for the matrices Gi, Gr 
in (2.11) and (2. 28), (2. 29), respectively. First, it is possible that F has more eigen- 
values inside the contour than the dimension of the space "H. For example, this 
is typical for characteristic equations of ordinary delay equations. However, this 
cannot occur with infinite dimensional spaces "H, IC which is our main concern here. 
Second, there may be a vector ^2=i ^nVn that is annihilated by all test functionals 
Wj or a functional PjWj that annihilates all test functions Vk (similarly for 

CSGEs). If the cigcnhmctions ?;„ and U7„ arc linearly independent and the data 
Wj,Vk are chosen at random we consider this case to be nongeneric. However, it 
is possible for nonlinear eigenvalue problems that eigenvectors belonging to differ- 
ent eigenvalues arc linearly dependent (sec [B, Sec. 4, 5] for such an example). The 
contour method can be extended to handle all these degenerate cases by evaluating 
higher order moments 

/ X''E(X)dX, p = 0,...2K -1. 

It is shown in [B, Lem.5.1] that it suffices to take K = J2n=il^";i under the 
conditions of Theorem 2.8. O 

2.4. Approximation of operators. We consider a sequence of approximate op- 
erators Fn G M{^1,F{Hn,I^n)),N G N, and study the errors for the linear system 
(2.3) and the spectrum <j{F) when F is replaced by F^. We will use the framework 
of discrete approximations (see [V76] and [V80] for an English reference) which 
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has the advantage that HnjICn are general Banach spaces (not necessarily sub- 
spaces of v., JC) connected to %, K only via a set of linear operators (not necessarily 
projections) 

PN-n^UN, qn-Ic^icn, N€N. 

In Section 4 we will apply the theory to boundary value problems on the infinite line 
when approximated by two-point boundary value problems on a bounded interval. 
In order to assist readers unfamiliar with the theory we impose conditions slightly 
stronger than necessary, and we try to avoid as many notions as possible from [V76] . 
Our assumptions are as follows: 

(Dl) There exist Banach spaces 'Hn,1Cn-,N e N and linear bounded mappings 
Pn e jC{'H,'HN),qN S £(/C, /Cat) with the property 

lim ||pArt;||-H„ = \\v\\n, v eU, lim ||gjvv||A:„ = e /C. 

JV— >-oo JV— >oo 

(D2) Given F e H(r2, J^{H, K,)) with p{F) ^ and Fn e H(f2, F{Un, JCn)), N e 

N with sup^gpj sup_;^g2 ]li<jv(A)|l < oo for every compact set C C O. 
(D3) -Fjv(A) converges regularly to F{X) for all A S in the following sense: 

(a) limjv^-oo \\Fn{X)pnv - qNF{X)v\\ =0, v eK, 

(b) for any subsequence vn € 'Hn,N G N' C N with \\vN\\nr,r,N G N' 
bounded and limN/gAr-^oo \\F]s[{X)vn — QnuWiCn = fo'^ some y € K,, 
there exists a subsequence N" C N' and avGH such that 
limN"9Ar->oo \\vn - PnvWun = 0- 

The term \\Fn{X)pnv — qNF{X)v\\ is called the consistency error since it measures 
the defect in the noncommuting diagram. 



(2.32) 



The following therorem summarizes convergence results from [V76, §3(3), §4(33)]. 

Theorem 2.10. Under the assumptions (D1)-(D3) the following assertions hold: 
( i) For any compact set C C p{F) and any v £ IC there exists an Nq e N such 
that the linear equation FN{X)yN = q^v has a unique solution t/jv = yjv(A) for 
N ^ No, X £ C, and the following estimate holds 

sup\\yNiX) - pNyiX)\\n^ ^ C sup\\FN{X)pNy{X) - qNFiX)yiX)\\, (2.33) 
xec xec 

where y{X) G H is the unique solution to F(X)y = v. 

(a) For any Aq G 'j(F) there exists Nq G N and a sequence Ajv G <j(Jjv), N ^ Nq, 
such that Ajv — Aq as N ^ oo. For any sequence Aat € (T(-FAr) with this convergence 
property, and associated eigenvectors G ^(i^jv(Ajv)), ll'^^ArllHiv = 1; ^^^^ following 
estimates hold: 



Xn-Xo\ < Cejv" (2-34) 
1 

"■0 

N ' 



inf ^\\v% - PnVo\\-Hn <C!e^, (2.35) 
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where ro — Hi is the maximum rank of eigenvectors that belong to Xq . The quantity 
En is a consistency error defined by 

En = max m&x\\FN(X)pNV — qNF(X)v\\, 

|A-Ao|^<5 veM 

where M. = spanjiij^j, : j = 0, . . . , /Zj, — l,p = 1, . . . , g} is the root space associated 
with Ao (see Definition 2.5) and 6 > is chosen sufficiently small. 

Remark 2.11. We note that our assumptions (D2) and (D3) imply regular conver- 
gence in the sense of [V76, §2(17)]. Moreover, it is easy to see that the convergence 
result in [V76, §3(3)] holds uniformly in A € C. Of course, for the latter result 
continuity of A F{X) is sufficient. O 

3. Relation to the Evans function 

In this section we relate the abstract approach from section 2 to the study of 
zeroes of the Evans function. Recall from the abstract setting formula (2.12), which 
we rewrite as 

E{X) = ^H,{X)+Ho{X), Xen\{Xi,...,X^}, (3.1) 

where we have introduced 

£,{X) = n7=i(A - A,), Hi{X) = Gi diag (n-=i,,^,(A - A,), j = 1, . . . , x) GJ. 

Note that ^^(A) may be viewed as an abstract version of the Evans function since 
its zeroes are exactly the eigenvalues of F{X). Moreover, Hq,Hi are holomorphic 
matrix functions such that Hi degenerates to a rank one matrix at every simple 
eigenvalue. 

In the following we make the above relation more explicit when the classical 
definition of the Evans function is used, that is as the determinant of a matrix with 
holomorphic columns determined from appropriate subspaces. In the; first step we 
set up a normalized Evans function that is independent of a perturbative situation. 

3.1. Definition of a normalized Evans function. In general the Evans function 
may be considered as a determinant that measures the coalescence of two families 
of subspaces depending holomorphically on a parameter. We define a normalized 
version of the function that is unique up to a sign. To begin, we introduce the class 
of matrices 

y^d,k = {p g l^d-k . dct(P^P) 1} (3.2) 

and recall the following elementary fact: Let U C be a subspace of dimension 
k. Then two matrices Pi,P2 e of rank k satisfy 7^(Pl) = Tl{P2) = Z7 if and 
only if there is an invertible matrix R G C'^''^ such that P2 = P\R. 

Definition 3.1. We call two matrices Pi,P2 G C''''^ equivalent and write Pi w P2 
provided P2 = PiR for some R G C'"''" such that det(i?) = -1 or det(i?) = 1. We 
will use notation [P] = {P\ G C'^'*' : Pi w P} for the equivalence class. 



Remark 3.2. If Pi, P2 G M'^'*' and 7^(Pl) = 7^(P2) then Pi w P2 holds automatically 
since 1 = det(P2"^P2) = det(P7Pi) det{Rf = det{Rf. O 
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Lemma 3.3. For any rank k projection 11 in there exist P, € C'''*' such that 

n = P*"^, $"^P = 7fe. (3.3) 

Moreover, if the inequality det(P"'"P) ^ holds for some P from (3.3), then it holds 
for any such P, and there exist Po, $o € C''''^ such that 

n = Po$^, PoeM'''^ $JPo = 4. (3.4) 

Proof. Using basis vectors in 7?.(n) and 7?.(n^) as columns, we obtain rank k ma- 
trices P,$ G C'^''= with -^(n) = 7^(P) and 7^(^^) = 7^($), respectively. Standard 
linear algebra shows A/'($^) = ??.($)-'- = 7^(11^)-'- = A/'(n), where means orthog- 
onal with respect to the duality pairing (sec (2.1) in Section 2). 

First wc obtain invcrtibility of <1>^P G C*''*': For c G C*^ the assumption <i>^Pc = 
implies Pc G 7^(P) nA/'($^)^ = 7^(^) n A/'(n) = {O} and hence c = since P 
has full rank. Second, we claim that 

n = P($'^P)-i$'^, P,$e (3.5) 

Indeed, the matrix = P($^P)^i$^ is a pro.joction with 7^(^) = 7^(P) = 7^(^) 
and 7V(n) = A/'($^) = .A/'(n), yielding H = ft and finishing the proof of (3.5). 
Letting 

$ = $(pT$)-\ (3.6) 

we arrive at the representation (3.3). 

If P, P are two matrices as in (3.3), then P = PR for some invcrtiblc R G C'^'*^ 
due to 7^(P) = n{P) = 7^(^). Now det(P^P) = det(P^P) det(P)2 proves the 
second assertion in the lemma. Finally, if det(P^P) ^ then we can replace P 
and $ in (3.5) by 

Po = Pdiag(det(P'^P)-i/2, 1, . . . , 1) e M''''= and $o = $(Po^5)"^ (3.7) 
Since Pq G Rf*'*, we arrive at (3.4). ■ 

Next consider two subspaccs of complementary dimension. We may write them 
as images of projections that both have the representation (3.3), or, perhaps, even 
the normalized representation (3.4). 

Definition 3.4. Let f7, F C he subspaces of dimension k and d—k, respectively. 
Pick any P G C'^''' of rank k such that ^{P) = U and Q £ C^'d-k ^^^^ ^ _ ^ 
such that 'R-{Q) = V . Then we call 

mniQ]) = {det(Piigi) : Pi G [P],Qi G [Q]} (3.8) 

the determinant set of the equivalent classes ([P], [Q]) from Definition 3.1. In 
addition, assume that P e M<^''=, Q e M'^''^-*'. Then we call 

V{U, V) = {det(P|Q) : P G M''''=, Q G M'^'''-'^, 7^(P) = U, n{Q) = V} (3.9) 

the determinant set of the subspaces U and V. 

Remark 3.5. Wc stress that the set X'([P]. [Q]) is defined with no additional as- 
sumptions on the subspaces U, V. Although this set does not depend on the choice 
of the representatives P, Q in the equivalence classes [P], [Q], it does depend on the 
choice of the equivalence classes. On the other hand, the set T>{U, V) is defined 
under an additional assumption on the subspaces U, V, but is uniquely determined 
by these subspaces. O 
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Our key observation is the following lemma which shows that the class (3.2) and 
the definitions (3.8), (3.9) have some significance. 

Lemma 3.6. Let U,V C be subspaces of dimension k and d — k, respectively. 
Pick any P € C''^ of rank k such that n{P) = U and Q e C'^''^-^ of rank d-k 
such that Tl{Q) = V. Then there exists z G C such that V^P], [Q]) = {z, —z}. The 
subspaces U and V are complementary if and only if z ^ Q. In addition, assume 
that there are Pq,Qo picked as above and such that det(Po'^Po) det((5Q Qo) 7^ 0. 
Then there exists z €<C such that V^U^V) = {z,—z}. 

Proof First take Pi e [P], Qi & [Q] with 7^(Pl) = U,n{Qi) = V and define 
z = det(Pi|<5i). Then take another pair P2, Q2 of this type and use Definition 3.1 
to write P2 = PiP, Q2 = QiS for some matrices R € C^'^ , S € C^-^^^-k g^^j^ ^YiaX 
det(P), det(5') e { — 1, !}• Therefore we conclude 

det(P2|g2) = det [{Pi\Qi) = ^det(P) det(5). 

This proves I'([P], [Q]) C {z,—z}. In fact wc have equality since —z is attained 
by taking R = diag(— 1, 1, . . . , 1), 5* = Id-k- The second assertion of the lemma 
is obvious. By normalizing Po,Qo as in (3.7), we may assume that Pq € M"^'*^, 
Qo e M'^'''~'=, and the third assertion holds by Remark 3.2. ■ 

Definition 3.7. A family of subspaces U{\) C C'', X E fl C C, is called holomor- 
phic in f2 if there exists a holomorphic family of projections 11 S M.{p., C'^''^) such 
that 7^(^(A)) = U{X) for all XgQ. 

Clearly, by the connectedness of fl the projections must be of constant rank and 
hence the subspaces are of constant dimension. 

The following lemma generalizes the representation (3.3) and the normalized 
representation (3.4) to holomorphic families. It is essential that the domain is 
simply connected. 

Lemma 3.8. Let be a simply connected domain in C. For any holomorphic family 
of projections 11 S 11(0, C''''') of rank k there exist functions P, $ G M{Q, C'^''^) such 
that 

n(A) = P(A)4>(A)^, $(A)^P(A) = 4, forallXen. (3.10) 

Moreover, for any P e H(n, €''■'=) such that 7^(P(A)) = TZ{U{X)) for all XeQ, the 
set 

A = {A e n : dct(P(A)^P(A)) = 0} (3.11) 
is discrete, may be em,pty, and is independent of the choice of P. Finally, there 
exist functions Po,^o G IH(f^ \ A, C''*^) such that 

n(A) = Po(A)$o(A)^, $o(A)^Po(A) = /fe, Po(A) e M''''^ for allXen\A. (3.12) 

Proof. By a result from [K, Sec.II.4.2] there exist functions P, $ G 11(1], C''''') such 
that 7^(^(A)) = Tl{P{X)), n{Il{Xy) = 7e($(A)) for all A e fi. This step uses 
simple connectedness. In the next step, as in the proof of Lemma 3.3, we normal- 
ize ^>(A) as in (3.6), and keep holomorphy. This proves (3.10). The set of zeros 
of the holomorphic function det(P(-)"^P(-)) is discrete. If Pi,P2 e H(n,C'*''=) are 
such that 7^(Pl(A)) = 7^(P2(A)) = 7^(^(A)) for all X € fl then there is a non- 
singular P(A) G C*^'*^ such that P2(A) = Pi(A)P(A); hence, det(Pi(-)^Pi(-)) and 
det(P2(-)^P2(-)) have the same zeros. If A e Q.\A then det(P(A)'^P(A)) ^ and we 
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normalize P{X) and $(A) as in (3.7), and keep holomorphy. Note that the square 
root has (up to a sign) a unique analytic continuation in any simply connected 
domain that does not contain 0. ■ 

The following theorem shows how one can assign to two families of holomorphic 
subspaces the usual Evans function which is equal to zero if and only if the subspaces 
are not complementary. This function is holomorphic on all of J7 see (3.13). It 
depends, up to a sign, on the equivalence classes from Definition 3.1 for the matrices 
formed by the basis vectors of the subspaces. In addition, one can assign to the 
two families of subspaces an Evans function that is unique, up to a sign, but is 
holomorphic on a smaller set O \ A(C/, V), see (3.15). We will call it the normalized 
Evans function. We refer to Remark 3.5 regarding the sets T>{[P{X)], [Q(A)]) and 
V{U{X),V{X)) used in (3.13), (3.15) below. 

Theorem 3.9. Let O G C fee a simply connected dom,ain and let U{X), V{X), X E fl, 
be two holom,orptiic families of subspaces of dimensions k and d — k, respectively. 
Pick any P G H(fi,C''''=) such that n{P{X)) = U(X), rank(P(A)) = k, and Q € 
H(0, C'^.'^-fc) such that 7^(Q(A)) = V{X), rank(Q(A)) = d-k for allXe Ct. Then 
the following assertions hold. 

(i) There exists a holomorphic function £ : — > C such that 

£{X) eV{[P{X)],[Q{X)]) for all A e O. (3.13) 

(m) £{Xo) = if and only if U{Xo) n y(Ao) 7^ {0}. 

{Hi) If £1 G H(f2, C) is any function satisfying (3.13) then either £1 = £ or 
£1 = -£. 

(iv) Moreover, the set 

A{U, y) = {A e n : det(P(A)'^P(A)) det(Q(A)"^Q(A)) = 0} (3.14) 

is discrete, may be empty, and is independent of the choice of P, Q picked as above, 
{v) Finally, there exists a holomorphic function £q : fl\ A{U, V) C such that 

£o{X) eV{U{X),V{X)) for all X&n\A{U,V), (3.15) 

and properties (i), {ii) hold for £q. 

Proof. By Lemma 3.8 we have a representation for the projections Ilc/jIIy e 
H(Q,C''''') associated with [/(A), ^(A) as follows 

nc/(A) = P(A)$(A)^, ny(A) = Q(A)*(A)^, A e O, (3.16) 

where P, $ e H(f2, C'*''^), Q, * e H(0, £d,d-k-^ and for all A e f2 

$^(A)P(A)=7fe, *^(A)Q(A) = 7d_fc. (3.17) 

Defining 

£{X) = dct(P(A)|Q(A)). Xen, (3.18) 

and using Lemma 3.6 proves assertions (i) and {ii). Now let £1 be any function 
with the properties of £. Then the quotient £i/£ is holomorphic on Cl \J\f{U,V) 
where the zero set is given by 

^^{U, y) = {A e O : V{[P{X)], [Q(A)]) = {O}} = {a e O : £iX) = 0}. (3.19) 

Since the quotient assumes only values ±1 there and Af{U, V) contains only isolated 
points it is a constant yielding (Hi). Assertion (iv) follows from Lemma 3.8, see 
(3.10). Also, by Lemma 3.8, see (3.12), we can choose functions Po,$o € W.{fl \ 
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A(i7,F), and Qo,*o G M{n\A{U,V),C'^-'^'-'') such that, in addition to (3.16), 
(3.17), the following normalization holds: 

Po(A) e M''''=, Qo(A) e M'^''*-'=, Xen\ A{U, V). (3.20) 

Defining 

£:o(A) =det(Po(A)|Qo(A)), Xen\A{U,V), (3.21) 

and using Lemma 3.6 proves assertion (v). ■ 

Remark 3.10. 1. An alternative proof of the theorem can be obtained directly 
from Lemma 3.6 without using Lemma 3.8 and [K, Sec. 11. 4. 2]. One first defines 
£, respectively. So locally via (3.18), respectively, (3.21) using a respective local 
normalization that is always possible (without simple connectedness). Then one 
shows that one can continue holomorphically along any curve in i7, respectively, 
O \ A{U, V) (since the set Af{U, V) defined in (3.19) is isolated) and the sign of the 
continuation is determined from Lemma 3.6. Since the domain is simply connected, 
the result follows from the Monodromy Theorem, see, e.g., [C, Cor. IX. 3. 9]. 

2. According to Theorem 3.9 the normalized Evans function is unique up to a sign 
on simply connected domains. The assumption of simple connectedness is crucial 
for the proof. The main ingredient which makes the Evans function unique up to 
a sign was to allow only matrices with det(P^P) = 1 in Definition 3.4. When 
homotopy arguments are applied to the Evans functions between different points 
in il, e.g. between and oo. it remains to be checked whether the values are from 
the same leaf that belongs to one of the signs. 

3. Normalization (3.7) shows that the normalized Evans function Sq can not be 
continued to A{U,V) since fo(A) ~ (A — Ao)~™/^ as A — >■ Aq at any point Aq € 
A(U, V) which is a zero of order m of the function det(P(-)'^P(-)) det{Q{-)'^Q{-)) 
holomorphic in i7. O 

In the next step we project vectors in onto their components in the subspaces 
U (A) and ^(A). Although this cannot work at the zeroes of the Evans function we 
insist that the main part of the projection stays holomorphic at the singularity. 
Recall the adjugate adj(A) of a square matrix A e C''"'^ given by 

(adj(A)),^ = (-1)^+^ det {A,,rny^-f:^^^l;:J > (3-22) 

(see e.g. [S, Sec. 4. 4]). It satisfies the identities ^adj(A) = adj(^)^ = dct{A)Ij and 
det(adj(A)) = det{A)'^~^. More importantly, by definition the adjugate preserves 
holomorphy, i.e. if A is in 11(0, C'^''') then so is adj{A). 

Theorem 3.11. Let the assumptions of Theorem 3.9 hold and let £ G M{^l,C), 
respectively, £q G H(0\ A([/, V^), C) he one of the two Evans functions, respectively, 
normalized Evans functions, determined there. Then there exist matrix valued func- 
tions yu,yv € IHI(0,C''''') with the following properties: 

£{X)Id^yuiX)+yv{X) for all A e 0, (3.23) 

n{yu{x)) = u{x), Af{yu{x)) = v{x), 

n{yv{x)) = v{x), Afiyv{x))^uix) ^w^u, [6.z^) 

n{yu{x))cu{x), jV{yu{x))DV{x), 

TliyviX)) CV{X), Af{yv{X))DU{X) ^w-u. K^.^o) 
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Conversely, the functions yu^yv G 11(17, C'^''') are uniquely determined by proper- 
ties (3.23) and (3.24). Similarly, for the normalized Evans function £q there exist 
yuW,yvW e B.{fl\A{U,V),C'^''^) with the same propeHies. 

Remark 3.12. Note that yu and yv behave almost like projections since (3.23), 
(3.25) imply 

yu{\)yu{\) = £Wyu{X), yv{\)yv{\) = £{\)yv{x), Aea (3.26) 

o 

Proof. For the projections Hu^^v S ]HI(r2, C'*''') associated with the subspaces 
f/(A),y(A), we take the representation (3.10) from Lemma 3.8 and partition the 
adjugate as follows 

adj (P(A)ig(A)) = (f [J|t) , A e O, (3.27) 

where R G H(r2, C''''^), S G H(17, C'''''~'=). With these settings wc define 

yu{\) ^ P{X)R{\y , yv{\) = Q{\)S{\y for (3.28) 
Conditions (3.23)- (3.25) hold for A G 17 as can be seen from the identities 

(P(A) Q{X)){^^^^-,^=£{X)I,= {J^^-,YpW W)- (3.29) 

Note that the first equality implies (3.23) while the second implies (3.25). In case 
£{X) ^ the matrices i?(A), S{X) are of full rank so that (3.24) follows. 

Now consider functions yu,yv S 11(17, C'*''') that satisfy (3.23), (3.24) in 17. 
If £{X) ^ then from the general form (3.5) and (3.24) we find representations 
%{\) = P{X)R{\Y, yv{X) = Q{X)S{Xy for some .R(A) e C'^'^5(A) e C<^''^-'=. 
Invoking (3.23) leads to 

(P(A) Q(A))(^||^J^)=£:(A)7,, 

and hence R{X),S{X) and i?(A),5(A) must agree due to (3.29). Thus % = yu, 
yv = yv holds in 17 \ N{U, V) and hence in all of 17 due to holomorphy. 

The last assertion in the theorem follows as above by taking the representation 
(3.12) instead of (3.10). ■ 

Next we consider the behavior of the matrices yu and yv near zeroes of the 
Evans function £. Let us first note, that the multiplicity of a value Aq as a zero of the 
Evans function corresponds exactly to the algebraic multiplicity of the eigenvalue 
Aq of the matrix pencil 3^(A) — (P(A)|Q(A)) from (3.18) defined in terms of root 
functions (see [MM, Prop. 1.8. 5]). For simplicity we consider only the behavior 
near simple eigenvalues and show that the matrices yu-,yv degenerate to rank one 
matrices. 

Theorem 3.13. Let the assumptions of Theorem 3.11 hold and let Aq G 17 &e 
a simple zero of the Evans function £{X) defined as in (3.18). Then there are 
nontrivial vectors vo,iuo such that 

^^{P{Xo)\Q{Xo)) = span{t;o}, ^((P(Ao)|Q(Ao))^) = span{«;o}. 
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Morover, with ^^o = ( ) we have the formulas 



yu{Xo) = £'{Xo)P{Xo)vo,iw^ = -S'{Xo)Q{Xo)vo,2W^ = -3^y(Ao) (3.30) 
and, if VqVq = 1, 

£'{Xo) = vj adj(P(Ao)|0(Ao))(P'(Ao)|Q'(Ao))t;o. (3.31) 

Finally, if X G il \ A(U, V) then the Evans function £ in assertions above can be 

replaced by the normalized Evans function £q . 

Proof From [MM, Prop. 1.8. 5] we have that the kernel of (P(Ao)|Q(Ao)) and of 
its transpose arc one-dinicnsional. Applying Keldysh's Theorem to the matrix 
pencil y{X) = (P(A)|Q(A)) shows that for some hi € H(Jli, C'^''*), Oi is some 
neighborhood of Aq, 

y{X)-^ = -T^vow^ + hi{X), Xefli. (3.32) 

A — Ao 

Here we have normalized wq such that 

1 = wJy{Xa)va = ujJP'{Xo)vo,i+w^Q'{Xo)vo,2- 
Comparing the singular parts of (3.32) and 



leads to 



't ] = voWq = ' J{ . (3.33) 



f^(Ao) \SiXo) ' ; " ~ \vo,2wl 
With this and = y{Xo)vo = P(Ao)wo,i + Q{Xo)vo,2 we arrive at 

yu{Xo) = P(Ao)P(Ao)^ = £:^(Ao)P(Ao)^;o,i«^o^ 

= -.?^(Ao)Q(Ao)^;o,2«^^ = -3^v(Ao). 

Finally, normalizing VqVq = 1 wc obtain (3.31) from (3.33) and the normaliz- 
ing condition. Note that (3.31) may also be derived directly by differentiating 
adj(3^(A))3^(A) — £o(A)/fc at A = Aq and then multiplying by Vq from the left and 
by Vq from the right. I 

3.2. Application to general first order differential operators. In this sec- 
tion we apply the previous results to general first order differential operators with 
matrices that depend holomorphically on the eigenvalue parameter. We consider a 
first order {d x d) matrix differential equation, 

y' = A{X,x)y, x€R, A e O, (3.34) 

and the respective pencil of first order differential operators 

F{X)y = ~y' + A{X,x)y, a; e M, A e fi, (3.35) 

and impose the following assumptions. 

Hypothesis 3.14. (i) Assume that the mapping 

n3 X^ A{X, •) - P(-) e L°^(R, C''^'') is holomorphic (3.36) 

for some matrix valued function B{ ) such that either B{-) G L^(M, C*''') or B(-) is 
bounded continuous with limx-).±oo B{x) = 0. 
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(ii) Assume that the differential equation (3.34) has for a/Z A G il exponential 
dichotomy on M+ with projections P+{X,x), x ^ 0, of rank k and on R_ with 
projections P-{X,x), x ^ 0, of the same rank k, i.e. for every X £ Q there exist 
C,a > so that the usual dichotomy estimates in (3.41) below hold. 

(Hi) Assume that (3.34) has exponential dichotomy on M for some X Gfl. 

Hypothesis 3.14(i) yields that the differential operators F{X) are defined on the 
A-independent domain, cf. [DL, Lem.2.8], [CL, Ch.3,4], given by 

n = {y€ L^R, &):ye ^Cioc(K, C''), -y' + B{-)y e L\R, &)}. (3.37) 

Moreover, F{X) is a bounded operator from H into the space 

)C = L'^{R,&), (3.38) 

when H is equipped with the graph norm 

\\y\\n = \\y\\h + \\-y' + B{-)y\\h- 

If B{-) is bounded, then H = iI^(M, C^), the Sobolev space. The completeness of 
H follows from Lemma A.l saying that H is embedded in the space of continuous 
functions vanishing at ±oo ,see also [DL, Lem.2.8]. 

We note in passing that the differential operator F{X), considered as an un- 
bounded operator in L^(M, C^) with the domain H, generates a strongly continuous 
semigroup {T*}t^o, called the evolution semigroup, see [CL], defined by the for- 
mula {T*y)(x) = S{x, X — t, X)y{x — t), a; G K, f ^ 0. Hero and below we denote by 
S'(a;,^, A), a;, ^ e M, the propagator (solution operator) of the differential equation 
(3.34). 

Remarks 3.15. (a) Hypothesis 3.14 (ii) holds if and only if the operator F{X) is 
Fredholm and its Fredholm index is equal to zero for all A G il. This is Palmer's 
Theorem, see [P84, P88], and also [BG], [Sa, Thm.3.2], [SaSl, Thm.2.6] and [LT, 
LP, SaSS] for its discussions and generalizations. In particular, the operator F{X) 
is invertible if and only if (3.34) has an exponential dichotomy on the entire line R. 

(b) Hypothesis 3.14 (iii) implies that p{F) ^ 0, cf. [BL, SaSS], [CL, Thm.3.17], 
and thus Hypothesis 3.14 yields Hypothesis 2.1 for pencil (3.35). 

(c) Since the coefficient of the differential equation (3.34) is holomorphic, the 
dichotomy projections P±{-,x), x £ M±, are holomorphic, see, e.g., [BL], [DL, 
Lem.Al], [SaSO, Thm.l] and further references therein. O 

Example 3.16. (Perturbations) A typical case where Hypotheses 3.14 are met 
occurs with operators of perturbation form 

[F{X)y]{x) = -y'{x) + {Ao{X, x) + B{x))y{x), a; e M, 

if the unperturbed operator Fo{X)y = —y'+Ao{X, x)y has an exponential dichotomy 
on M for all A € and with projections of rank k. Then the assumption B G 
L^(M, C'^''') guarantees that the exponential dichotomies on half lines hold for the 
perturbed operator (see, e.g., [Co, Prop. 4.1], [GLM, Lem.2.13]). This case applies 
to the Schrodinger equation in Section 5 below. O 

A more specific situation occurs when one linearizes a one-dimensional PDE 
about a traveling front and re- writes the respective eigenvalue problem as the first 
order ODE system (3.34). Then its coefficient can be assumed to stabilize at ±oo. 
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Example 3.17. (Traveling Fronts) Consider the case of picccwisc constant matrices 
Apc{X,x) = A+{X) for x^O and Apc{X,x) = A_{X) for a; < 0, (3.39) 

where the matrices A± (A) satisfy the following properties: 

(a) ^±(A) analytically depend on A G fJ, 

(b) ^±(A) have no purely imaginary eigenvalues, 

(c) rankP4^(^) = rankP4_(;^) for the Riesz projections Pa±{\) corresponding to 
the part of the spectrum of A±{X) located in the left half plane. 

Under these conditions the (unperturbed) differential equation y' = Apc{X,x)y 
has exponential dichotomy on M.± and is Fredholm of index zero by Palmer's result 
cited in Remark 3.15(a) above. This corresponds to fact that A does not belong 
to the essential spectrum of the underlying differential operator that appears when 
one linearizes the PDE about the front, cf. [Sa]. Again, as in the previous ex- 
ample the perturbed differential equation y' = (^pc(A,a;) + B{x))y inherits the 
dichotomies on half lines if either B S L^(R, C'^''') or B{-) is bounded continuous 
with limj;^±(x) B{x) = 0. In both cases we are in the system class described in 
Hypothesis 3.14. We cite [BL, LT, SaS8] for further references. O 

As described in Section 2, we consider an inhomogenous equation 

F{X)y{X) = v G L"^. (3.40) 

As above, let S{x,^,X),x,^ € M, A € 17 denote the solution operator of F{X) (the 
propagator of the differential equation (3.34)), and let P±{X,x) be the dichotomy 
projections for (3.34) on ]R_|_,]R_ from Hypothesis 3.14(ii). Then we have for some 
C, a > 0, 

5(x,^,A)P±(A,e)-P±(A,x)5(x,e,A) for x,^eR±, 

||^(x,e,A)P±(A,e)i| s;Ce-«(^-«), for ^^x,x,^eR±, (3.41) 

||5(a;,^,A)(/-P±(A,e))KCe-"(f-^), for x^^,x,^gR±. 

Let £ be the Evans function from Theorem 3.9 (i) with respect to the equiva- 
lence classes [P(A)], [Q(A)] for any choice of holomorphic functions P,Q such that 
7e(P(A)) =7^(P+(A,0)) and7e(0(A)) = A/'(P_(A, 0)), i.e. 

£(A) e2?([P(A)],[0(A)]), Xen. (3.42) 

Alternatively, let So be the normalized Evans function from Theorem 3.9 {v) with 
respect to the subspaces U{X) = 'R{P+{X,0)) and V{X) = J\f{P-{X,0)) in C^, i.e. 

£o{X) e V{U{X), V{X)), XGfl\ A{U, V). (3.43) 

We solve (3.40) for A e p{F) in a standard way by using Green's operators 



{g+iX)v){x)= G+{x,^,X)md^, x^O, 
Jo 

(g. {X)v) {x) = [ G_ {x, X)v{Od^, x^O, 

J —OO 



(3.44) 



acting on functions : M — >■ C'', with the kernels given by 

ri^.y^-j S{x,^,x)P+{x,0 o«;e^x, 

U+(x,l;,A)-<^ ^(a;,e,A)(P+(A,0-/) 0^x<t 

/ 5(a;,C,A)P_(A,0 ^^x^^O, 
U-[x,^,A)-<^ 5(a;,C,A)(P_(A,0-/) xK^^O. 



(3.45) 
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Due to the exponential dichotomies the operators Q+ , Q~ have uniform bounds in 
all spaces Lp(R,C''), 1 < p ^ oo, see, e.g., [Co] or [CL, Sec.4.2]. The following 
piecewise defined function gives the general solution of (3.40) on both half lines 
(we write v+ = v\R_f_, V- = ■!)|r_ for short), 

^^ = [^(.^,0, A)r,+ (A) + ((?+(A)i,+) (x), x > 0, 
^'^ \-5(a;,0,A)r7_(A) + (g_(A)i)_)(x), a; < 0, ^' ' 

where ?7+(A) e U{X) and »7-(A) e V{X) are arbitrary. The function defined in (3.46) 

is a solution of the differential equation (3.34) if the left and right limits of j/(A) at 
zero coincide. Note that if y± G ACioc(K±, C'') then the function y defined by 

y+{x), X > 0, 
y-{x), X <0, 

is in ACioci^, if ^^'^ on^Y if V+i^) = V-i^)- Therefore we have to find vectors 
r/+(A) e U{X) and 77_(A) e ^(A) such that 

7j+{\) + r;_(A) = (g_(A)z)_)(0) - {g+{\)v+)iO) =: [0]o. 

By Theorem 3.11 the sought for vectors arc given by 

1+^^^ = SlXj^u^^'^^^^o, V-W = -^yviXMo- (3.47) 

Inserting this into (3.46) gives the solution formula 

,.A.^-/ MX)S{x,0,X)yu{X)[i]o + {g+{X)i+){x), x>0, 

The same formulas hold with £ replaced by Sq provided X £ Q \ A{U,V). It is 
convenient to introduce the operators g{X) acting on functions t; : M ^ C'' by 

{g{X)v){x) = { (3.49) 
V V / ;w {g-{X)viu_){x), x<0, 

and the matrix valued function G by 

G{X,x) = l ^i^'Vm;^m ^ n M 

' ^ \ -s{x,o,x)yv{x), x<o, ^ ^ 

so that if ?;o e is a given vector then G(A, ■)vo : K — >■ C''. 

Using these notions in (3.48) and inserting thcini into (2.4) finally leads to the 
expression (3.51) below for Ejk{X),j = 1, . . . ,m, k = 1, . . . ,i. To formulate the 
result, recall that the operators F{X) act from H into /C , see (3.35), (3. 37), (3. 38). 
Also, note that Ti C L^(M, C^) C n' such that a function w e L^{R, C'') defines on 
H a linear functional by {v\v) = {w,v)r; here and below for w,v £ L^(R, C*) we 
denote (i«,w)r = J^^w{x)^ v{x) dx. Following Section 2, we now choose linearly 
independent functions Vk G L^(M, C*), k = !,...,£, and linearly independent func- 
tions Wj e Z/^(IR, C*), j = 1, . . . , m, viewed as elements of H'. Thus, the discussion 
above can be recorded as follows. 

Theorem 3.18. Assume Hypotheses 3.14, and let all eigenvalues of the operator 
pencil F inside fig be simple. Then the matrix E{X) G C"*'^ from the contour 
method (2.4), (2.5) satisfies the following formula 

Ejk{X) = -^{wj,G{Xr)[vk]o)R + {wj,G{X)vk)R, A en, (3.51) 
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where the vector [£'fe]o = {G-{X)vk){Q-) - (^+(A)f;fe)(0+) E denotes a jump 
quantity at x = 0, and the operators Q, Q± and the matrix valued function G are 
defined in (3.44), (3.49) and (3.50). For \ & Vl\K{U,V) formula (3.51) holds with 
£ replaced by the normalized Evans function So ■ 

Formula (3.51) shows how to express the abstract terms in (3.1) through integral 
kernels and the Evans function (or the normalized Evans function) for first order 
systems. 

To conclude this subsection, we summarize our results for the operator pencil 
(3.35). In particular, we apply Theorem 3.13 to recover the singular part of the 
function Ejk{ ) near a simple eigenvalue A„, n = 1, . . . , x. Recall that -P(A), (3(A) € 
M{n, C''*^) are chosen as in (3.16) with Uu{X) = P+(A, 0) and ny(A) = I-P- (A, 0), 
and the Evans function is defined by £{X) = det (P(A)|Q(A)). 

Theorem 3.19. Assum,e Hypotheses 3.14. The following assertions are equivalent. 

(i) Ao is a sim,ple eigenvalue of the operator pencil F (3.35); 

(ii) Ao is a simple root of the Evans function £; 
(Hi) dimAA(P(Ao)|Q(Ao)) = 1; 

(iv) dim(7^(P+(Ao,0))^A^(P_(Ao,0))) =1; 

(v) There exists a unique up to a scalar multiple exponentially decaying at ±oo 
solution V of (3.34); 

(vi) There exists a unique up to a scalar multiple exponentially decaying at ±oo 
solution w of the adjoint to (3.34) equation {z^)' = ~z^A(X,x)^. 

Moreover, ifv{0) denotes the initial value of the exponentially decaying solution in 
{v), then v{0) e 'R{P+{Xq,0)) n.A/'(P_(Ao,0)) if and only if 

v{0) = P(Ao)wo,i = -Q{Xo)vo,2 (3.52) 

for a vector vo = (^^"'^^ from Af{P{XQ)\Q{XQ)) . 

In addition, assume that all eigenvalues A„, n — 1,...,k, in are simple, 
let Vn,wJ^ he the solutions described in assertions (u), [vi) above for each A„, and 
normalized as indicated in (2.8), let {vkYk=i' {'^jlj^i linearly independent func- 
tions in L^(]R, C^) chosen as indicated in Section 2. Then the singular part of the 
function Ejk{-) defined in (2.4) in the framework of the abstract Keldysh theorem 
(2.10) is given by the formula 

^ 1 

^jrW = Z]"^ T-kwhVnWkwn^vuV. |A- A„| < 1, n= l,...,x. (3.53) 

n=l ^ 

Proof. The three assertions, £:(Ao) = 0, 7^(P+(Ao,0)) n 7V(P_ (Aq, 0)) 7^ {0}, and 
A/^(P(Ao)) 7^ {0}, are equivalent by Theorem 3.9 and by the dichotomy assump- 
tions in Hypothesis 3.14. Moreover, 7^(P+(Ao,0)) n 7V(P- (Ao, 0)) and 7V(P(Ao)) 
are isomorphic via the map v(fS) 1— > w(-) = 5(-, 0, Ao)w(O). Thus to see the equiv- 
alence of the first four items in the theorem it suffices to show that the sub- 
spaces 7^(P+(Ao,0)) n A/'(P-(Ao, 0)) and A/'(P(Ao)|(5(Ao)) are isomorphic as in- 
dicated in (3.52). Let ^;o e M{P{Xq)\Q{Xq)). Since 7^(P(Ao)) = 7^(P+(Ao,0)) and 
7e(Q(Ao)) = AA(P_(Ao,0)), we have P(Ao)wo,i e 7e(P+(Ao,0)) and Q(Ao)«o,2 e 
A/'(P- (Ao, 0)), and by the choice of uo we have P(Ao)t'o,i = — Q(A„)z;o,2. Then 
v{^) from (3.52) belongs to 7^(P+(Ao,0)) n A/'(P- (Ao, 0)). Conversely, if v(fS) € 
7e(P+(Ao,0)) n j\^(P-(Ao,0)) then v{0) = P{Xo)vo,i and v{0) = -Q{Xo)vo,2 for 
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some vo i € C'', wo 2 e Letting vq = ( '^°'M yields vq e 7V(P(Ao)|0(Ao)), as 

required. To begin the proof of (vi), we remark that dimA/'(P(Ao) |(5(Ao)) is equal 
to dimA''((P(Ao)|Q(Ao))^); also, the following identities hold: 

M{{P{Xo)\QiXo)V) - 7^((P(Ao)|0(Ao)))^ = n{P{Xo))^ n 7^(Q(Ao))^ 

= 7^(P+ (Ao, 0))^ n Af{P. (Ao, 0))^ = Af(P+{Xo, 0)^) n 7^(P_(Ao, 0)^) 
= n{i - P+(Ao, 0)^) n Af{I - P_(Ao, 0)^). 

Since (3.34) has the exponential dichotomy P±(Ao,0) if and only if the adjoint 
equation has the exponential dichotomy / — P±(Ao,0)^ (see, e.g., [BG, Lem.4.5], 
[DL, Lem.2.4], [Sa, Rem.3.4]), it follows that the subspaces j\^((P(Ao)|(5(Ao))"^) and 
J\f{F{Xo)^) are isomorphic via the map w{0) 1— )• w{-)'^ , where w is the exponentially 
decaying solution of the adjoint equation. 

It remains to show (3.53). For each n = 1, . . . , x let Vn{0) be the vector from 
(3.52), let Vn{x) = S'(.T, 0, A„)i;„(0) be the decaying on both and IR_ solution 
of (3.34), and let w„(0) span 7V((P(A„)|Q(A„))^). We obtain from (3.30) that 

yu{Xn) = f'(A„)t;„(0)«;„(0)^ and yv{Xn) = -f '(A„)t;„(0)ti;„(0)^, 

and from (3.50) that G(A„,a;) = 5'(A„)t;„(a;)w„(0)^. Using (3.51) and £{X) = 
S'{Xn){X — Xn) + 0{\X — A„p) yields the singular part 

A - A„ 

Note that 



..{OVlvkjo =wj{0) (^y° S{0,^,X)P.{X,Oi>k{Od^ 

S{0,^,X){P+{X,0-I)vk{Od^ 



Next we observe that 



(0)3(0, ^,x)p.ix,o, e^o, 

(0)5(0, C,A)(/-P+(A,e)), OO 



solves the adjoint equation of (3.34), is continuous at 0, and decays exponentially 
in both directions. Using this in (3.54) finally leads us back to the singular part 
determined by the abstract Keldysh theorem in (2.10), where the term {wnjVk) is 
understood as the integral {wn,Vk)R- B 



4. Convergence of eigenvalues for finite boundary value problems 

In this section we provide error estimates of the eigenvalues obtained by the 
contour method in Section 2 when the boundary value problems (3.40) are solved 
approximately on a bounded interval. In the first step we analyze the error of 
the boundary value solutions themselves, and in the second step we discuss the 
implications for the contour method. 
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4.1. Estimates of boundary value solutions. Using the setting and notation 
from Section 3.2 we consider the all-line boundary value problem 

F{X)y{X) = -y'{X, •) + A{X, .)y{X, ■) =v G L\R) 

for various values of A and v. Our main assumption is the following. 

Hypothesis 4.1. Let Hypotheses 3.14 hold and assume that the dichotom,y ex- 
ponent a > in (3.41) is uniform for all X & ^l and that the dichotomy projec- 
tions P±{X,x) given in Hypothesis 3.14 (ii) are asymptotically constant, that is, 
limcc-).±oo -P±(A,a;) = P±(A). 

We approximate (3.40) by a sequence of boundary value problems on finite in- 
tervals Jn = [x^,x^], N E N, with —x^, oo as ^ oo: 

where R-,R+€ H(Q,C<*''^) are given matrix valued functions. We allow boundary 
conditions that are nonlinear in the eigenvalue parameter in order to cover so-called 
projection boundary conditions which lead to fast convergence of solutions of (4.1) 
as —x^,x^ — >■ 00, sec (4.16) below. 

We will show how (4.1) fits into the framework of Section 2.4 and apply Theorem 
2.10 to obtain error estimates. Our approach is largely based on [BR] where the 
case of smooth coefficients and domF(A) = i?^(IR, C*), the Sobolev space, was 
analyzed. For any interval J C R introduce the Banach space 

Hj = {y e L^J,^) : y e AC\oc{J,C''),-y' + By e L^J,C'')} (4.2) 

with norm \\y\\lij = ||y||i2(j) + \\ - u' + By\\'l2^jy Note that Hr agrees with H 
from (3.37). Using Lemma A. 2 it is easy to see that the space Hj is complete in 
this norm (for later reference wc collect further properties oiHj in the Appendix). 

The spaces T-Ln and /Cat from Section 2.4 are defined by "Hjv = "Hj^ and /Cat = 
L'^{Jn,C'^) X C<^ with norms \\y\\nN and ||(w,r)|||^ = + \r\'^. The spaces n, 

K are taken as in (3.37), (3.38) and mapped into %n, K^n by 

PNy = y\jM and qNV = {v\j„,Q). (4.3) 

Obviously these mappings are linear, bounded uniformly in A'^ and satisiy condition 
(Dl) from Section 2.4. 

By Hypothesis 3.14 (i),(ii) wc have F e m{Vl,F{U,K.)) since F(A) is Fredholm 
for each A in and the map O 9 A H> A{X, •) - B{-) e L°°(M, C'^''') is holomorphic. 
Moreover, p{F) ^ by Hypothesis 3.14(iii), and (D2) follows from the next lemma. 

Lemma 4.2. Under the above assumptions on A,R-,R+ the operators F^ are in 
H(f2, J^('Hjv,^Ar)) and sup j^^f^ sup xec ll-^Jv(A)|| < oo for every compact setC C ft. 

Proof. Let y gTIn and A € O. Then, by (3.36) and Lemma A. 2, 

\\FN{x)y\\l, = \\-y' + A{x, Ml^j.) + \R-Wyi^-) + R+Wy{^+)\' 

< "^{Mn, + ll(^(A, •) - B(-))2/lli^(7.)) + 2{\\R-{X)f + \\R+iX)f) Ibllioo 

< c{\\A{x,-) - B{.)\\U + WR-Wf + \\R+m')\\yrn.- 

From the holomorphy of A,R^,R-^- wc obtain uniform bounds for ||Fjv(A)|| on 
compact sets C C as well as holomorphy of A Jjv(A) for all N £ N. Finally, 
the Fredholm property is a well-known fact for finite boundary value problems. ■ 
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For the application of Theorem 2.10 it remains to verify (D3). Let Vf.{X) be a 
basis of the range of -P-(A) and let V^(A) be a basis of the kernel of f+(A). 

Proposition 4.3. Let Hypothesis 4.1 hold and assume that matrices R±{X) for all 
A e satisfy 

det ((i?_(A)Kf (A)|i?+(A)V7(A))) ^ 0. (4.4) 
Then -Fjv(A) converges regularly to F{X) for all A e O. 

Proof. Let 2/ e A e O. Then 

\\F^{X)pNy-qNF{X)y\\l, < 2{\\R.iX)f + \\R+{X)f){\y{x'!)f + \yK)f), (4-5) 

where the right-hand side converges to as A/' — )• oo by Lemma A.l. This proves 
part (a) of (D3). 

Let A e O be fixed. Consider a subsequence j/at G Hn, N e N', with bounded 
||yiv||-Hjv assume there is v € IC with limAr->.oo \\FN{X)yN — qNv\\icN = 0- Set 
{vN,rN) '■= Fjq{X)yN G K-N and note that j/at can be written similarly to (3.46): 



VNix) 



G+{x, 0, A)yAr(O) - G+{x, x^,X)yN{x^) + [Gf (A)«Ar](a;), x ^ 0, 



-G-{x,0,X)yN{0) + G-{x,x^,X)yN{x^) + [G^ {X)vn]{x), x < 0, 

(4.6) 

where G± are defined in (3.45) and, cf. (3.44), 



{G!^iX)vN){x)= G+{x,^,X)vN{Od^, 
Jo 

{G'!{X)vn){x)= [ G-{x,^,X)vN{Od^, 



xl^x^O, 



a;^ < a; < 0. 



By Lemma A. 2 the boundedness of ||2/]v||-Hn implies boundedness and thus com- 
pactness of the sequence {yN{^))New C C^, i.e. limN"3jv->-cx) J/iv(0) = yo for some 
subsequence N" C N' and some yo € C*. Let 

/^■^ ^\y+{x) = G+{x,Q,X)yQ + [G+{X)v\-fi^]{x), x ^ 0, 
\y-(x) = -G-(x,0,A)yo + [e-(A)«|Kj(x), x<0. 

Now we follow verbatim the proof of [BR, Thm.2.1] until [BR, (2.11)] to conclude 

Note that this step is crucial. It uses the determinant condition (4.4) as well 
as the representation (4.6) and the exponential dichotomies. It remains to prove 
\\yN —PNy^^UN — >■ as N" 9 A'' — >■ oo, for which the arguments in [BR] do no longer 
apply. 

By construction y± G ACioc(K±, C'') n L'^{R±,C'^) and 

-y'± + A{X, ■)y± = W|M^ (4.9) 
holds in L'^{R±,C'^). Therefore, the function defined by 

\-y'_+B{x)y-{x), x<0, 
is in L^(M, C^) and satisfies z = v + (-B(-) — A{X, ■))y. Using this we obtain 

l|-2/jv+5(-)2/JV-^|jjvlU=(Jiv) = \\-y'N+B{-)yN-{v+{B{-)-A{X,-))y)\j^\\L^j^) 
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^ \\FN{X)yN - QNvW^, + \\A{X, •) - B{-)\\L'^\\yN - y\jJmj.), (4.11) 

where the right-hand side converges to zero as N" 9 A/' — >• oo by (4.8) and our 
assumption. Without loss of generality we may assume ./a? D [—1,1] for all e N". 
Repeating the estimate (4.11) with [—1,1] instead of Jjv shows that (t/jvl[-i,i])7VGN" 
is a Cauchy sequence in Therefore, its limit, which coincides with 

is an element of ^C([— 1,1]). This allows us to conclude z = —y' + B{-)y from 
(4.10), so that equations (4.8) and (4.11) prove our final assertion. ■ 

The above results show that the abstract convergence result, Theorem 2.10, 
applies: 

Theorem 4.4. Let the assumptions of Proposition 4.3 hold. Then for any compact 
set C C p{F) n and any v e K. there is Nq & such that for all N ^ Nq and 
A e C the boundary value problem (4.1) has a unique solution yN{X,-) £ Hn- 
Furthermore, for some constant C, independent of v, 

snpWyNiXr) - PNy{X,-)\\n. ^ C snp\R-{X)y{X,x^) + R+{X)y{X,x^)\, (4.12) 
Aec Aec 

where y{X, ■) CzH solves (3.40). 

From the well known decay \y{X, x^)\ — >■ as iV — >• oo (e.g. see the proof of [BL, 
Thm.3.2]) estimate (4.12) implies that the solutions yjv(A, •) of the finite interval 
problems converge uniformly in A G C to the solution of the problem on the line. 

We will now concentrate on the differential equation (3.34) with the coefficient of 
the special perturbative structure which appears in the case of traveling fronts with 
asymptotic hyperbolic rest states, see Example 3.17. Specifically, let us consider a 
first order operator of the form (3.35) where 

A{X,x) = ^pc(A,;r) + B{x), x e R, (4.13) 

with Apc{X,x) defined in (3.39). We impose the following assumptions. 

Hypothesis 4.5. The differential equation (3.34) withA{X,x) from (4.13) satisfies 
Hypothesis 3.14 and Hypothesis 4.1 with the uniform exponential estimate, 

||P±(A,ar)-P±(A)|Kce-«l^l, x€R±, (4.14) 

for all A e fi, where a is the exponent from (3.41) of the exponential dichotomy on 
M-t for (3.34). The projections P±{X) depend analytically on A € O. 

A typical situation where Hypothesis 4.5 is satisfied, occurs when the matrix- 
valued function A{X,-) = Apc{X,-) + B{-) is continuous, assumptions (a) - (c) in 
Example 3.17 hold, and there is c > such that for all A e fi, 

\\A{X,x)-A_{X)\\ < ce-"l^l, a; < 0, \\A{X,x) - A+{X)\\ < ce""'^!, a; > 0, (4.15) 

where a is the exponent of exponential dichotomy on M± for the constant coefficient 
equations y' = A±{X)y. Then the differential equation (3.34) with A{X,x) as in 
(4.13) has exponential dichotomy on 1R± for all A e O and the roughness theorem 
[BL, Thm.A.3] implies (4.14). Thus Hypothesis 4.5 is satisfied provided (3.34) has 
an exponential dichotomy on M for at least one A G f2. 

Under Hypothesis 4.5, projection boundary conditions in (4.1) are a suitable 
choice, because they always satisfy (4.4) by construction. For convenience, we recall 
the definition of the projection boundary conditions, see [B90] for more details. 
Since the limits P±{X) depend analytically on A, there are analytic bases V^{X), 
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respectively, V^{X} of 7^(P±(A)), respectively, A/'(P± (A)) (e.g. [K, Sec.II.1.4]). Let 
us split the inverse matrix composed as follows: 

The projection boundary conditions are then given by the boundary matrices 

By construction {R_{X)V1{X)\R+{X}V^{\)) = Id, i.e. (4.4) is satisfied, and 

i?_(A)(J-P_(A))=0, i?+(A)P+(A)=0. (4.17) 

We now use Theorem 4.4 to establish our main convergence result for the matrix 
in (2.4) and the integrals in (2.5) when the underlying boundary value problems 
are solved on a finite interval. We denote by = A^g(]R, C*) the set of finite, 
compactly supported, valued Radon measures on R. By Riesz's Theorem, e.g. 
[F, Thm.7.17], and Lemma A.l, C H' for H from (3.37). U w e H' is given 
by /U e Ml, we write {w,v) = J^v''^ {x)dfi for v e H- We approximate iu = fi G 
Ml C "H' on a finite interval J by its trace w\j = defined through {w\j,v) = 
JjV^{x)diJ, for all v G Hj. Obviously, 

{w\j,Pjv) — {w,v) = / v''^ {x)dfi = 0, if J D supp(/i) and w e 'H. (4-18) 

Jm\j 

Example 4.6. Two standard examples for w G H' , given as /i G M^: 

(1) If /X = CiSxo for some Xq & ^, i & {l,...,d}, then {w,v) — j {x)dii = 
v{xq)^ Ci = t'i(.To) is the j'th component of v evaluated at xq. 

(2) If fjL has density / G L\^^{^, C*) with respect to the Lebesgue measure, then 

{w,v) = ^^v^{x)f{x)dx. O 

After these preliminaries we define and estimate approximations E^{X) of the 
function E{X) from (2.4) by solving finite interval boundary value problems. 

Theorem 4.7. Let F from (3.35) satisfy Hypothesis 4.1, and let condition (4.4) 

hold for boundary matrices R±{X). Moreover, let T C p{F) H be a contour and 
assume linearly independent elements wj G H' , j = 1, . . . ,m, defined by fxj G M^, 
and linearly independent functions Vk & K,, k = 1, . . . ,t which are hounded and have 
compact support. Then the finite interval approximation E^{X) of E{X), defined 
by (cf (4.1); 

FN{X)v^{X) = {vk\j^,Q),k = !,...,£, 

^^(A),fc = (A)) = / (A)^dM,-, j = 1, . . . ,m, 

J Jn 

satisfies 

sup ||£;^(A) - £;(A)|| < ce-"™^-^-'^^>, (4.19) 
Aer 

with a being the dichotomy exponent in from Hypothesis 4.1 and c a uniform 
constant. 
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// F additionally satisfies Hypothesis 4.5 and Fn is given with the projection 
boundary conditions defined via (4.16), then (4.19) improves to 

sup ll^^(A) - ^(A)|| < ce-2"™^-^-'0. (4.20) 

As a corollary we obtain estimates for the approximate matrices (cf. (2.5)) 

Df = ^.J^X^E''{X)dX, j = 0,l. (4.21) 

Corollary 4.8. Under the assumptions of Theorem 4.7 the following estimates 
hold 

\\Do - D^W ^ ce-"™{-^->0, \\Di - D^W ^ ce^^'^''^^-^- (4.22) 

In the case of projection boundary conditions the constant a improves to 2a. 

Proof of Theorem 4.7. Throughout the proof, c is a generic constant. For all A e 
r C p{F) the differential equation (3.34) has an exponential dichotomy on M with 
a uniform exponent a and projections P(A, x) = P±(A, x) that depend holomorphi- 
cally on A and satisfy limx-).±oo -P(A, x) = -P±(A), see [BL, Thm. A. 5]. With these 
projections the Green's function reads 



S{x,tX)P{X,0, x^^, 
J{x,^,X){P{X,0-I), x<^, 

and the solution yk{X) of F{X)yk{X) = Vk is given by (see [BL, Thm.A.l]) 



G{x,i,X) 



ykiX,x) = / G(.T,e,A)^fe(OdC. 
Jr 

There is A^o G N with Jn D supp(/Uj) for all N ^ Nq and j = 1, . . . ,m. Using (4.3) 
and (4.18) we find for A^ ^ A^o 

\{wj,ykiX)) - (%|.7„,yf (A))| 

< (wj, yfe (A)) -(%!,/„, pAryfe (A)) + {wj\j^,pNyk{X) - yk {X)) 

^ c pNykiX) - yj^ {X) . (4.23) 

Since Vk has compact support and is bounded there is a constant c such that 
\vk{0\ < ce-2"lf for all ^ G R and A; = 1, . . . This is used to bound the right 
hand side of (4.12): 

|i?_(A)j/fc(A,x^)Kc||i?_(A)|| / ^ e-"^^--^%{0\d^ 

J —OQ 

POO 

+ c\\R.{X){P{X,x'l)-I)\\ / e-«(€---)|%(0|rfC 
<c r"e-«l-^-«le-^«l€ld^ (4.24) 

J — OO 

poo 

+ c||i?_(A)(P(A,a;^!)-/)||e«^- / e-««e-2«l«l 
< ce'^"^- + c\\R-{X){P{X,x^) - I)\\e"^- . 



28 



W.-J. BEYN, Y. LATUSHKIN AND J. ROTTMANN-MATTHES 



A similar estimate holds for \R+{X)yk{X, x^)\. Since the estimates are uniform in 
A e r, we obtain 

sup|i?_(A)t/fe(A,a;^) + i?+(A)yfc(A,0| < ce-«™{---'-+>, (4.25) 
Aer 

by Theorem 4.4 this proves (4.19). 

If Hypothesis 4.5 holds, the projections P{X, x) of the exponential dichotomy on 
the whole real line can be chosen to satisfy (4.14), again sec [BL. Thm. A. 5]. For 
projection boundary conditions we then find from (4.14) and (4.17) 

\\R.{XmX,x^) - = \\R.{XmX,x^) - P_(A))|| < ce--. 
Summarizing, we can bound the right hand side of (4.12) as follows 

sup\R_iX)ykiX,x^) + B,+ {X)yf,iX,x^)\ < ce'^" '<>, 
Aer 

which gives the desired improved order of convergence. I 

Remark 4.9. It is not difficult to weaken the assumption of compact support for 
Vk and Wj. For example, the proof in (4.24) shows that it is sufficient to have 
\vk{0\ < ce-2«ISI for ^ e M and fc = 1, . . . O 

4.2. Estimates of eigenvalues. We analyze the error of the numerical algorithm 
from Section 2.2 when the matrices Dq, Di are replaced by their approximations 
£'0', Di satisiying the estimates (4.22). 

As before, let Dq be of rank x and let — VoSoVFg, be the short form of 
its singular value decomposition (SVD), cf. (2.16). In the following we consider a 
small perturbation Dq G C™'^ of Dq with the full SVD 



£'0= (1 



Vo L ° ~° . (4-26) 




where Vo e C™'^,^! e C™'"-^, M^o € C^'", Wi € C^'^"^, and Sq € C^'^ contains 
the X largest singular values of Do- Instead of computing the eigenvalues of 

D = Vo*DiWo^o^ (4.27) 

in (2.18), we use (4.26) and compute the eigenvalues of 

D = Vo*D,WoY:o\ (4.28) 

where Di is a small perturbation of Di . The following lemma shows that the eigen- 
values of D approximate those of D with the order of the original perturbations. 
In order to apply the perturbation theory from [St] we use the Frobenius norm 
\\D\\p = tv{D*D), the spectral norm ||£>||2 and the Hausdorff distance 

distij(Mi, M2) = max( sup inf |2 — y|, sup inf |^; — y|), Mi,M2CC. 
zeMi yeMa ^gMa y&Mi 

Lemma 4.10. Let Do,Di € C"*'^ be given such that rank£>o = >^ o-nd such that 
the matrix D from (4.27) has only simple eigenvalues. Then there exist So > and 
C > such that the spectra of D and D from (4.28) satisfy 

distH {(T{.D),a{b)) ^C{\\Do-Do\\f + \\D^-D4f), (4.29) 
provided \\Do — Do\\f + \\Di — Di\\f ^ ^o- 
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Proof. In the following C denotes a generic constant depending on eo but not on 
Do, and s := || A) - A)||f + \\Di - A||f < £0- 

Let us extend Yq = Va and Xq = Wa to unitary matrices Y = (Yq Fi) e C"'"* 
and X = {Xo Xi) e C^^^ and introduce E e C'"^'^ such that 

So Oj^^_j^ \ 77. / -£'11 -£'12 



Y=(yo Yi)=Y 

X = (xo Xi) = X 



(4.31) 



\'Jm->c,>c 'Jm->c,e->c/ \^21 -C'22/ 

with conformal partitioning. We require 4eo < o'min = iniiij=i,...,j<: cTj and obtain 
7 := (||£i2|||. + ||^2i||f)'^' < ||-E||f < £ < £0, 

2£ < 2eo < < CTmin - V^EO < CTmin - ||-Ell||2 - II-E22II2 =: S. 

Therefore [St, Thm. 6.4] applies and yields Q e C""^'^, P G C'"^'^ with 

{\\Q\\% + \\P\\%Y^' ^^^^^e, (4.30) 

SO that the unitary matrices 

7 -Q*\ f{I + Q*Q)-'/^ 0^,m-. 

,Q I ){ 0™_^,^ (7 + QQ*)-l/2 

/ -P*\ /(/ + P*P)-V2 o^,^_^ 

^ A (/+pp*)-v2 

transform Dq into block diagonal form 

Y*DoX = ^" G C™'^. 

The matrices An and A22 can be written as 

All = (/ + Q*g)'/' (So + Ey^ + i;i2P) (/ + p*p)-^/2 , 

A22 = (/ + QQ*)'/' (£^22 - £;2iP*) (-T + PP*)-^/' . 
Prom the estimate (4.30) we find 

(/+g*Q)±l/2_ J ^^j|Q||2^^£2 
2 

and a similar estimate holds for (/ -f QQ*)'^'^!'^. Using (4.31), (4.32) this leads to 
\\% - V^Wf ^ Ce, ||lo - W^p < Ce, |lSo - Anib ^ Ce, HAaajb Ce. (4.33) 

By [GK, II Cor. 2.3] and the last two inequalities in (4.33) the sets of singular 
values sing All of An and sing A22 of A22 satisfy 

min(sing All) > Cmin - Ce, max(sing A22) < Ce. 

Decreasing eq further we find that sing An and singA22 are disjoint and hence 
sing All contains the k largest singular values of Dq. 

Therefore, Ti{XQ) is the invariant subspace of DqDq corresponding to the x 
largest eigenvalues of DqDq and coincides with TZ{Wo). Similarly, TZ{Yo) = TZ{Vo). 

Then the matrices To = X^Wq G C^"" and So = Y^Vq G <C'^''^ are unitary and 
satisfy XqTo = Wq and Fo'^o = Vq- Moreover, Eq = S^Yq*DoXoTo = S^AnTo, so 
that the matrices 

D = Vo*D,Wof:o^ = S^Yo*D,XoToT^A^,^So = S^Y^* DiXqA^^' So 



(4.32) 
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and 

D = Y*D^XoA^^ 
are similar and have the same spectrum. Now estimates (4.33) imply 

\\b- D\\f ^Ce. 

Since simple eigenvalues depend analytically on the matrix, we finally obtain for 

some C > 

distjj {u{D),a{D)^ = distjj {u{D),a{D)^ < Ce. (4.34) 
This finishes the proof. ■ 

Combining this result with Corollary 4.8 shows that using boundary value prob- 
lems for the computation of the point spectrum is a robust method and leads to 
exponential convergence with respect to the length of intervals. 

Theorem 4.11. Let the assumptions of Theorem 4.7 hold and let all eigenvalues 
of F inside the contour T he simple. Then there is C > such that the following 
holds for all intervals [x^,x^] with min{—x^,x^} sufficiently large. The set 
of the eigenvalues of the approximate pencil from (4.1), computed by the method 
from Section 2.2 using the approximations Dq , from (4.21) instead of Dq,Di, 
satisfies the estimate 

disin (cr^, (y{D)) ^ Ce-"'"'"^-^-'^?>. (4.35) 
In case of projection boundary conditions the constant a im,proves to 2a. 

Remark 4.12. We note that the simplicity of eigenvalues was only used in the very 
last step (4.34) of the proof of Lemma 4.10. Similar to Remark 2.9 convergence of 
spectra as £ — )■ still follows in the general case from the perturbation theory in 
[K]. But now the rate is e~i^ in (4.34) where /i is the maximal algebraic multiplicity 
of eigenvalues inside V. Correspondingly, the rate a in (4.35) deteriorates to ^. 

5. The Schrodinger operator on the line 

We consider the eigenvalue problem for the one dimensional Schrodinger opera- 
tor, H, 

{H -X)u = 0, H = -d^/dx"^ + V{x), a; e K. (5.1) 

Here, the real valued potential V satisfies V <E L"'^(K), the domain of H is given by 

dom{H) = {u& L^{R) :u,u' e ACioc(M), -u" + VuG L^(M)}, 

and we assume that A e f2 = C \ [0, oo). Since V e -£'^(IR), the essential spectrum 
of H is equal to [0, oo), and the discrete spectrum consists of no more than finitely 
many negative simple eigenvalues > Ai > A2 > • • • > A,,, sec, e.g.. [RSIV, 
Sec. XIII. 3], [CS, Sec. XVII. 1.3]. The eigenvalue problem (5.1) can be written as the 
first order differential equation 

y' = A{X, x)y, A{X, x) = A{X, 00) + B{x), xgR; (5.2) 

here and below we denote 

In Subsection 5.1 we consider the linear operator pencil of second order differential 
operators F^^^\X) = H—XI, see (5.1), acting from the space H = dom{H) equipped 
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with the graph norm into the space /C — L^(M, C^). In Subsection 5.2 we consider 
the nonlinear operator pencil of first order differential operators F^^' (A) = —dx + 
A{\, •), sec (5.2), (5.3), acting from the space H as defined in (3.37) into the space 
IC — L^(M, C^). Our objective is to illustrate the construction of the matrix E from 
(2.4), (3.51), (3.53) and also its computation via some approximation arguments 
related to the boundary value problems on finite intervals, that is, to the equation 

-u"{x) + V{x)u{x) - \u{x) = 0, .T e [.T^. (5.4) 

equipped with appropriate boundary conditions at the endpoints x^, x^ satisfying 
— oo and a;^ +oo as N ^ oo, and to the boundary value problems 
for the first order differential equation (5.2). In the current section we do not 
assume exponential decay of the perturbation. Although we offer some explicit 
formulas for the matrix E and its approximation E^ in terms of certain solutions 
of the differential equations (5.1), (5.2), (5.3), we emphasize that they are mainly 
of theoretical value as our general approach in practical applications is not to use 
these formulas but instead to construct E^ by solving boundary value problems 
on finite intervals numerically. 

5.1. Second order differential operators. We consider the linear operator pen- 
cil F(^^'(A) = H — XI with H as in (5.1). Our main tool will be the Jost solutions 
u±{X,x), X & R, X G Q = C \ [0,oo), of the second order Schrodinger differen- 
tial equation (5.1) which are uniquely determined as the solutions of the Volterra 
integral equations 

/•±<X 

u±{X,x) = e^'^ X-^/^sm{X^/^{x-^))V{^)u±{X,Od^,x€R. (5.5) 

Jo 

Here and everywhere below we choose the branch of the square root such that 

Im(A-'^/^) > for A G f7, in particular, e*^ ^ ^ as .t +oo. It is well known 
that the Jost solutions satisfy the asymptotic boundary conditions 

lim e'F*^'^'"=?i±(A,.T) = 1, (5.6) 

a;— ^±oo 

they are holomorphic functions of A G f2, for A < they are real valued and positive 
for ±x sufficiently large, see, e.g., [CS, Chap. XVII]. The Wronskian 

W(A) = W{u-,u+) = U-{X,x)u'_^_{X,x) - u'_{X,x)u+{X,x), a; e M, A e O, (5.7) 

of the Jost solutions is equal to zero precisely at the points A„ e fi, the isolated 

eigenvalues of the Schrodinger operator H, where the exponentially decaying at 
+00 solution M_|_(A„, •) is proportional to the exponentially decaying at — oo solution 
u_(A„, •) with a nonzero constant c„, that is, when 

U+{Xn,x) = CnU-{Xn,x), X € R, Cn e C \ {0}, n = 1, . . . , X. (5.8) 

We refer to [W87] for the general theory of Sturm-Liouville differential opera- 
tors (see also [W05] for a brief but exceptionally readable account). In particular, 
due to y e Z/^(]R) the Schrodinger operator H is in the limit point case at ±00, 
the Jost solutions u± are L^-solutions at ±00, and thus the resolvent operator 

(f(")(A))~^ = {H - XI)-^ for A e ^1 \ {Ai, . . . , A;^} is the integral operator with 
the kernel 

1_^ U+(A,x)._(A,a -oo<C^.<+oo, 
^' W{u+,u_)\u-{X,x)u+{X,$,), -00 < a; < ^ < -Fog. ^ ' 
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Therefore, if Wj,Vk <E L^(K) are chosen as indicated in Section 2, that is, such that 
rank ((%(•), u±(A„,-))r)"„^i > x, rank ((u±(A„, •),%(•))«) ^^^i > x, 

then the matrix ^(A) = ((«;,(•), ((J^^^H^)) "'^^fe) (•))R),"^ii from (2.4) is given by 
^jfcW = ^;777 \ Wj{x)u+{\,x) u-{X,^)vk{^)d^dx 



yv{u+,U-j J-oo Jx 



(5.10) 



Integrating E{X) from (5.10) over the contour F from Section 2.2 we thus obtain 

X = rank_Do and formulas (2.18), (2.19) for the cigcnvahics of H. 
We now equip equation (5.4) with self-adjoint boundary conditions 

u(a;^)cosw± — u'(a;^) sina;± = 0, with some u)± e [0,7r), (5-11) 

and define the operator in L'^{[x^, x^]) by = -cf/dx'^ + V{x) with 

dom(iJ^) = {u eL^{[x^,x^]) -.u^u' G ACioc{[x^ , x^]) , -u" + Vu& ^^([a;^, x^]) 
and both boundary conditions (5.11) hold}, 

cf. [W05, Sec. 7]. Let u±{X,x) denote the A^-dependent solutions of the Schrodinger 
equation (5.1) each of them satisfying one of the respective initial conditions 

M±(A, xl) = e*^'''"£ sin w±, u'±{X, x^) = e'^"^^l cos w±. (5.12) 

Since m+(A,x), respectively, m_(A,.x) satisfies the boundary condition (5.11) at x'^, 
respectively, x^ we conclude (see, e.g., [W05. p. 84]) that the resolvent operator 
(f ("'^)(A))"^ = {H^ - Xiy^ is the integral operator with the kernel 

T,N,y^t:^_ 1 /M+(A,a;)M_(A,0, X^i^^^X^X^, 

W{u+,u.)\u_{X,x)u+{X,0, x^^x<^^x'l. ^ ' 

In the following we consider the restriction pnu = u\^^n^^n^ (cf. (4.3)) as an 
operator from ^^(K) into L"^ {[x^ , x^]) . Then it is known from [W05, Thm. 7.1] 
that {F^'''^\X))-'^pn converges strongly to {F'^"^{X))-^ in L'^{R). This is called 
generalized strong resolvent convergence of to H in [W05] (here, one imbeds 
L'^{[x^,x^]) into L^(R) by setting functions equal to zero in R \ [x^,x^]). It 
follows that the matrix E{X) from (2.4) for the operator pencil F'^^^)(A) can be 
written as the limit of the matrices E^{X) defined via the approximative operator 
pencils F("'^)(A). 

Proposition 5.1. Assume V G L^{M.) and let E{X) be defined as in (5.10). Then 

N 

' w^{x)({F^"-^\X))-^v^){x)dx) ' , 
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and we denote Wj' — pNWj, = pNVk- Similarly to (5.10), using (5.13) the 
matrix E^{X) can be computed by the formula 



W{U+,U-) J^N J^N 

N N 

J-— wf{x)u.{\x) u+{X,Ov^iOd^dx. 



(5.14) 



Wc note in passing that hniTv^oo W(m+, m_) = C{X)W{u^, m_) where the factor 
C(A) can be exphcitly computed and is equal to zero precisely at the eigenvalues 
of the operators and defined as H^. = —d'^/dx^ on L^((— oo, a;^]) and 
L'^{[x^ , oo)) with the domain determined by the respective boundary condition in 
(5.11), see [LS, Thm.3.3]. 

5.2. First order differential operators. W ith A{X,x) as in (5. 2), (5. 3) we con- 
sider the operator pencil F^^\X)y = —y' + A{X,x)y and follow step-by-stcp the 
constructions in Subsections 3.1 and 3.2 culminating in formulas (3.51) and (3.53). 

First, we need to choose the projections on the subspaces of the initial values 
of the solutions of (5.2) exponentially decaying at +00 and —00 and construct 
their representation (3.16) and the respective Evans function (3.18). The matrix 
A{X, 00) for A e has no pure imaginary eigenvalues and the differential equation 
y' = A{X, oo)y has the exponential dichotomy on R with the dichotomy projection 
being the spectral projection P{X, 00) := P+(A) = P-W for A{X, 00) corresponding 
to the eigenvalue iA^/^. We recall that if A G then the eigenvalues ±iA^/^ of 
A{X,oo) satisfy Rc(iX^/^) < < Re{-iX^/^). 

Since the perturbation B in (5.2) satisfies ||S(-)|| S L-'^(M) by the general theory 
in e.g. [BL, Co, GLM] there exist dichotomy projections 

P±(A,x) ^ S{x,0,X)P±{X,0)S{0,x,X) (5.15) 

on IR.± for the perturbed equation y' = A{X, x)y such that the dichotomy subspaces 
'R,{P+{X,x)) = span{t/+(A, a;)} and J\f{P-{X,x)) = span{?/_(A, a;)} are uniquely 
determined while their direct complements j'V(P+(A, a;)) and TZ{P-{X,x)) are arbi- 
trary. Here and below we use notation 

y±(A, a;) = fej^' x€R, X€fl, (5.16) 



^u'±(A,a;) 

for the (2 x 1) vector solutions of (5. 2), (5. 3) which correspond to the Jost solutions 
u±{X,x) of (5.1) defined via the Voltcrra equations (5.5). Also, given a vector 

V = (^^^ € denote = ^ and = {v'^)'^ = (— ^2 vi) so that for 

any two vectors v,w € we have det(t'|'u;) = v^w = —v'^w'^; thus, if w,w e 
are linearly independent then the projection on span{w} parallel to span{u>} is the 
matrix {w-^v)~^ 



vw 



± 



To make a choice of the dichotomy projections P±{X,x) one is tempted to let 

P+(A,0) = (W(«_,^i+))-i2/+(A,0)j/_(A,0)^, 

/-P_(A,0) = (W(«+,u_))-i2/_(A,0)y+(A,0)^ 

and then use (5.15) to define P±(A, x) for x G R±. This choice, however, is not sat- 
isfactory as these projections are meromorphic in ft with the poles precisely at the 
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eigenvalues A„ while the constructions in Subsection 3 require holoniorphy. Another 
choice is to normalize u± by letting u±{\,x) = (mj.(A, 0) + w^(A, 0))~^/^u±(A, a;) 
and then replace u± and y± in (5.17) by u± and y±. This choice yields holomorphy 
of P±(A, x) mn\A, where A = {A : (uHX, 0) +u^(A, 0)){u'i{X, 0) +w'^(A, 0)) = 0}, 
and admits the normalization of the type (3.4), (3.12) leading to the construction of 
the normalized Evans function Sq. However, the great disadvantage of this choice 
of P±{X,x) is that these dichotomy projections are not asymptotic as x ^ ±oo to 
the spectral projections oi A{X,oo). 

We now construct the holomorphic dichotomy projections that are asymptotic 
to the spectral projections at infinity. Our main tool will be the solutions j/^(A, •), 
y+{X, •) on R+ and yt{X, •), y" (A, •) on R_ of the differential equation (5.2) satis- 
fying the asymptotic boundary conditions 

lim e-*^'^'^y1_(A,a;) = V, lim e'^''"^y" (A, a;) = w, (5.18) 

X— >-+oo X— >-+oo 

lim e*^'^'^?/f.(A,a;) = w, lim e-'^''^''yl{X,x) = v, (5.19) 

where v = ^j^^i/s^' w = ^ iX^/^ eigenvectors of A{X,oo) — (^^^^ q 

such that A{X, oo)v = iA^/^v, A{X, oo)w = — iA^/^w. The spectral projections 
of the matrix A(X,oo) are given by the formulas P(A,oo) — (2iA^/'^)^^vw^ and 
/ - P(A, oo) = -(2iAi/2)-iwv-L. 

Since ||-B(-)I1 G -^"'^(IK-)i the existence of the solutions yJ"(A, •) on satisfying 
(5.18) and the solutions yl'"(A, •) on K_ satisfying (5.19) is guaranteed by the 
celebrated Levinson theorem from asymptotic theory of differential equations (see, 
e.g., [CoL, Probl.III.29] or [E, Thms.1.3.1, 1.8.1] and also [GLM, Rem.7.11]). These 
solutions are obtained by solving certain inhomogeneous Fredholm type integral 
equations on semilines with holomorphic integral kernels and inhomogeneities. In 
fact, inspecting the proof of the theorem, cf. [E, Sec. 1.4] or [GLM, Thm.8.3], we 
observe that the solutions y'±."{X, •) are holomorphic in A € il and limiting relations 
(5.18), (5.19) hold uniformly in A on compact subsets of fl. 

We will now discuss how the solutions y^"(A, •), yl'"(A, •) satisfying (5.18), (5.19) 
are related to the solutions ?y±(A, •) defined in (5.16) via the Jost solutions u±{X, •) 
of (5.5). In fact, the solution y'^{X, •), respectively, 2/l(A, •) is uniquely determined 
by the first limiting relation in (5.18), respectively. (5.19) since due to (5.6) it is 
precisely the solution y+{X, •), respectively, 2/-(A, •) defined in (5.16): 

y'±{X,x) =y±{X,x), xeM.±, X e fl. (5.20) 

The solution 2/!|i:(A, •), respectively, 2/"(A, •) is not unique and can be changed by 
adding a summand proportional to y+(A, •), respectively, y_(A, •). A convenient 
choice of the solutions y±{X, •) is furnished by the formulas 

yl{X,x) = — :ry-{X,x) - 2^ — —y+{X,x), X e M+, 

W(u-,u+) ^c„(A-A„) 

(5.21) 

= W(u-,u+) "+^^'^^ - (A^A„) ^ ^ 

where y±(A, •) are defined in (5.16), A„, n = 1, . . . , k, are the eigenvalues of H, 
the constants c„ are taken from (5.8), and we denote by p„ the residue at A„ of 
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the function 1/W{X) for the Wronskian W{X) defined in (5.7). As we will see in 
a moment, the solutions j/j_(A, •), y'j.{X,-) defined in (5. 20), (5. 21) satisfy (5.18), 
(5.19). Using these solutions at a; = we let 

P+(A,0) = (2zy/2)-iy;(A,0)?y^(A,0)^, 

/ - P_ (A, 0) = -(2zAi/2)"iyl (A, 0)2/!l (A, 0)^ 

and then use (5.15) to define P±{X,x) for x € M.±. By a direct calculation we also 

have P_(A,0) = (2iAi/2)-iy« (A, 0)yi(A, 0)^. 

Lemma 5.2. Assume that V e L"'^(M) and A G = C\[0, oo). Then the projections 
P+{X,x) = {2iX^/^)-^y'^{\,x)yl(X,x)^, x&R+, (5.23) 
P^(X,x) ^{2iX^/^)-^y'^(X,x)yl(X,x)^, x G M_, (5.24) 

defined via formulas (5.15), (5.22) are holomorphic in Q and satisfy 

lim P±(A,x) = P(A,oo). (5.25) 

X— >±oo 

Proof Formulas (5.18), (5.19) for the solutions j/j_(A, •), y'!^{X,-) defined in (5.20), 
(5.21) follow from (5.6) and the relations limj,_j.±oo e^^^^^^^u'j_{X,x) = ±iA^/^ and 

lim e±*^'^'^u^(A,a;) = W(A)/(2iAi/2), lim e^'^'^'^u' (A,a;) = tW(A)/2, 

(see, e.g., [LS, Lem.3.1]). Also, computing the residues of the RHS of (5.21) and 
using (5.8) we see that the solutions yj_(A, •), y^{X,-) defined via (5.20), (5.21) are 
holomorphic in A G Jl. In the remaining part of the proof we concentrate on the 
case of M+ as the arguments for ]R„ are similar. The Wronskian of the solutions 
y^(A, •) and t/+(A, •) is x-independent and by (5.18), (5.19) we infer that 

yliX,x)^yl{X,x) = {e'^''SliX,x))^{e-'^''S^+iX,x)) ^ w^v = 2iX'^^ 

as x ^ +00. Therefore, 

yliX, x)-^yX{X, x) = 2iA^/^ for all x e M+ (5.26) 

and thus (5.23) is a holomorphic projection. But y^{X,x) — S{x,0, X)yf^_{X,0) by 
the definition of 5'(x, 0, A) and y^{X,x)-'- = y^{X,0)-'-S{x,0, X)~^ since w{x) = 
y+{X, 0)-'"5'(x, 0, A)~^ satisfies the adjoint equation w' = —wA{X, x)^ and therefore 
should be of the form w = y-^ for a solution y of the equation y' = A{X, x)y. Thus, 
(5.23) is in concert with (5.15) and (5.22). Using (5.18), (5.19) again we have 

y;(A,x)2/«(A,.T)^ = {e-'>^"'\%{\x)){e'>^"''^yl{\x))^ ^ vw^ 
as X ^ +0O yielding (5.25). ■ 

We are ready to identify the ingredients in the representation (3.16) of the pro- 
jections n[/(A) = P+(A,0) and ny(A) = I — P_(A,0) on the holomorphic families 
of one dimensional subspaces U{X) = span{y+(A, 0)} and V{X) = span{2/_(A, 0)}. 
Indeed, (5.20) and formulas (5.23), (5.24) show that (3.16) holds with 

P(A) =y;(A,0), $(A) = (2iAV2)-i(j;«(A,0))t, (5.27) 

Q(A) = -2/1 (A, 0), *(A) = (2iAV2)-i(y«(A,0))t, (5.28) 
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where the normaHzation (3.17) has been shown in the proof of Lemma 5.2, see 
(5.26). Therefore, due to (5.20) the Evans function (3.18) is given by 

£{X) =det{y+{X,0)\y-{X,0)) =W{u-,u+), (5.29) 

while the matrices yu{X),yv{X) in Theorem 3.11 are given by 

yu{X) = y+iX, 0)j/_(A, 0)^, yviX) = 2/_(A, 0)y+(A, 0)^. (5.30) 

Furthermore, the Green's kernels (3.45) of the Green's operators (3.44) are given 
as follows: 

r (n^^ \\ = ^_ /y;(A,a;)2/!^(A,0^, O^^^x, 
2iXy^\yliX,x)yl{X,0^, 0^x<t 

r( .^\\_ ^ iy^{X,x)yt{X,0^, ^^x^O, 
"^-^^^^'^f 2iAV2 |y^(A,ar)2;«(A,0^, x<^^0, 

while the function (3.50) and the vector [v]o = {g^{X)v){0-) - {g+{X)v){0+) from 
Theorem 3.18 are computed as follows (recall that y±{X, •) = y±{X, •) by (5.20)): 

G(Xx) = [^+^^'''^^-^^'^^^' 

[y_(A,x)y+(A,0)^, x < 0, 

[v]o = f yl{X,0)y.{X,0^v{Od^- yli>^,0)y+{>^,O^md^)- 

To give a compact formula for the singular part of the RHS of (3.51) it is convenient 
to introduce, for given functions WjV G Z/^(M, C^), the following (2 x 1), (2 x 1) and 
(2 X 2) matrices: 

y+{X, O^wi^) / y+{X, ^^vi^) ' 

, V{X) =1 I , (5.31) 




II^ y- (A, e)^5)(0 dcj V/°^ y- (A, O^m 

1 /-2/_(A,0)^|/!^(A,0) 2/_(A,0)^2/!l(A,0) 



2iAV2 V-y+(A,0)^y!^(A,0) y+(A,0)^y!i(A,0) 



(5.32) 



Plugging (5.21) in (5.32), a short calculation reveals the following concretization of 
formulas (3.51) and (3.53) for the case of the Schrodinger operator. 

Theorem 5.3. Assume V G L^(R) and use linearly independent functions Wj,Vk G 
iy^(M, C^), j = 1, . . . ,m, k = 1, . . . ,i, in (5.31). Then the singular part of the RHS 
of (3.51) can be expressed as follows: 

^ {wj,G{x,-)[vk]o)R = Wj{xfy{X)Vk{X),xen. (5.33) 



£{X) 



As in (3.53), using formula (5.33) the singular part of Ejk{X) near the eigenvalue 
X„, n = 1, . . . , X, of the operator H can be computed as follows: 

oo oo 

£;-^(A) = ^^( J y_{Xr^,^)-^wj{Od^y[ J y+{Xr^,^)^VkiOd^). (5.34) 



EIGENVALUES AND BOUNDARY VALUE PROBLEMS 



37 



Wc recall that at A = A„ the solutions y+(A„,-) and y_(A„,-) from (5.16) are 
proportional, see (5.8), and that the residue pn of the function 1/W(A) at the point 
A„ is given by 

_ 1 /■ dX 

- 2^ 4 wMvwvjy ^ ^ ^ 

for a sufficiently small circle 7„ centered at A„, sec (5.7). 

We now consider approximation (5.4) of the Schrodinger equation on the finite 
segment [x^,x^]. Let u±. denote the Jost solutions corresponding to the truncated 
potential defined on R by V^{x) = V{x), x e [x^,x^] and V^{x) = 
otherwise, that is, the solutions of the truncated Volterra equations 

/■^± 

u^{X,x)=e^'^ X-^/^sm{X^/^{x-^))V{^)u^{X,^)d^,x€R. (5.36) 

Jo 

Similarly to (5.16), we denote by 

y±(^'^)={(ui)^{x%)^ ^eR, xen, (5.37) 

the corresponding solutions of the first order diff^erential equation (5.2), (5.3). It is 
easy to see from (5.36) that these solutions satisfy the following conditions: 

y^(A, ) = e^^'''< V, y^(A, x^) = e"*^'''-- w, (5.38) 

where v,w are the eigenvectors of A{X,oo) defined after equations (5.18), (5.19). 
We will now identify the boundary conditions as required in (4.1) and (4.16). 

Since A+{X) ^ A_(A) = A{X,oo), we have P+(A) = P_(A) = {2iX^/^y^vMV'^ , 
TZ{P+{X)) = span{v}, A/'(P+(A) = span{w}, and thus the discussion in Section 4 
leading to (4.16) yields 

^+(^^=^Gav^ -i)' J)- ^^•39) 

In particular, due to (5.38) the solutions y±{X, •) from (5.37) satisfy the boundary 
conditions 

R-{X)y{X,x^) + R+{X)y{X,x^) = (5.40) 

used in (4.1) to define the approximate operator pencil Fjp{X). 

Wc are ready to formulate the final convergence result of this section for the 
Schrodinger case. Let j/^(A, •) from (5.37) be the solutions of the equation (5.2), 
(5.3) satisfying (5.36), (5.38), and choose linearly independent functions Wj,Vk € 
L^(IR, C^), j = 1, . . . , TO, k = !,...,£. Similarly to (5.34), at the eigenvalues A„, 
n = 1, . . . , >c, oi the operator H we consider 

oo oo 

</"'W= j ( j yt{>^n,O^Mi)d^), (5.41) 

— oo —CO 

where is defined by (5.35) with the Jost solutions ^^^(A, •) replaced by u± (A, •). 

Theorem 5.4. Assume V € L-'^(R). Then the singular part (3.53) of the matrix 
E{X) associated with the operator pencil F^^(X) on the whole line is the limit as 
N ^ oo of the matrices £J^'^'"S(A) in (5.41) associated with the operator pencil 
F^{X) on [x^,x^] with the boundary conditions (5.40). 
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Proof. Since the solutions y+(A, •) and 2/+(A, •) from (5.18) arc linearly independent 
solutions of (5.2), (5.3) and j/+(A, •) is a solution of the same equation, we may 
choose constants a^, such that y!^{X,x) = a^y'^{\,x) + /3^y'^{\,x) for all 
X e Solving the last equation for a^, (3^ at x = x^ and using (5.38), (5.26) 
yields 

^ _ det{y^{\,x^)\yl{\x^)) _ det{^\e'^'''<yl{\x^)) 

-2av^.5 .A. _ det(e-^^^^^<y;(A,0|v) 

P+ - _2ai/2 ^ 

as, N CO due to (5.18) (we recall that Re(iA^/^) < 0). A similar argument shows 
that if z/^(A,a;) = a^j/i(A, x) + /3^y!^(A, x) then 

a^^l, e^^^'''"-/?^ ^ asA^-^oo. (5.43) 

Computing Wronskians and using (5.20), (5.42), (5.43) we conclude that 

>V(y^ (A, ■),y^ (A, •)) = (A, •) + y!^(A, •), 2/;(A, •) + P^yliX, +)) 

converges to >V(j/-(A, •), y+(A, •)) as TV — )■ oo (and even uniformly in A on compacta 
in f2 because the convergence in (5.18), (5.19) is uniform on compacta). It follows 
that pn as N ^ oo. 

It remains to show that the integral terms in (5.41) converge to the respective 
integral terms in (5.34). The latter fact follows from the assertions 

\\y±i\-)-y±{\-)\\mR±n[x«.xf])^0 asiV^oo (5.44) 

using the Cauchy-Schwartz inequality. To prove (5.44) for M_)_ (the argument for 
M_ is similar), we recall that y+(A, •) = j/+(A, •) by (5.20), and then use the repre- 
sentation 2/+ (A, •) = a^?/+(A, •) + (A, •) to estimate 

lb+(A,-) -2/+(A,-)llz,2([o,^jv]) < |l-Q!+| ||2/+(A,-)||z,2([o,^iV]) 

+ \Pn\{ P e-^^^'^"'"^^ ■ e^^'^'^"''^^\yl{X,x)\Ux)'^^ 
Jo 

< |1 - |b+(A,.)|U=(M,) +c|/3;v|(e-^^^('^^'^K - l)'^' 



since the function e*'^^^^^?/+(A, cc) is bounded on R+ due to (5.18). Now (5.42) 
implies (5.44) finishing the proof of the theorem. ■ 

6. Numerical experiments for the FitzHugh-Nagumo equation 

We apply the contour-method to investigate spectral stability of traveling waves 
in the FitzHugh-Nagumo system (FHN). A traveling wave in FHN is a solution of 

Q = u" + CU' +U - - V, 

(6.1) 

= cv + $(u + a- bv). 

For the standard parameter values a = 0.7, b = 0.8, $ = 0.08 (sec [M],[BL, § 5]) 
one finds both, a stable pulse with speed c ~ —0.812 and an unstable pulse with 
speed c w —0.514. We choose the latter one and denote it by {u,v). 
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Linearization about this pulse leads to the linear operator 



u\ { u" + cu' + u — u^u ■ 

V 



cv' + ^u — ^bv 



For spectral stability one has to analyze the location of the spectrum of this oper- 
ator. The eigenvalue problem reads 

{XI -C) (:) = (''^ mL^m,C'). 



V J 

From the dispersion relation one can show (e.g. [BL, § 5] or [RM, § 7]) that there 

is no essential spectrum in = {ReA > —0.064}. Thus, we may use the circle 
r = {A S C : |A — 1| = 1.05} for the contour method. Equation (2.3) becomes 

(A7 - £) (^^l^j) = Vk in L^{R, C^). (6.2) 

For the computation we choose % as a random linear combination from the 2M- 
dimensional space 

Km = span{ (^^^^ , : j = 0, . . . , M - l} c L\R, C)^ 

lOj 

ifj{x) = max{0, 1 - \x-Xj\}, Xj = -5 + ^_ ^ , j = 0, .. .,M- 1. 

More precisely, Vk = J2jL~o^ {^lifj + Cfe(0, </5j )"""), with ^1 and independent 
and normally distributed random variables. 

With z = {yl,{yly jVk)^ we rewrite (6.2) as a first order system (see Section 
3.2), 

Xzi - CZ2 -zi+ u^zi + ^3 I + I ~vl I in L"(R, C^). (6.3) 
i$zi + i(A + $6)z3 

For the approximation (4.1) of (6.3) on a bounded interval we choose projection 

and periodic boundary conditions, which both satisfy (4.4). For these boundary 
conditions, we showed in Theorem 4.7 exponential error estimates (with a twice 
as good rate for the projection boundary conditions), when the Vk are compactly 
supported, a property shared by the linear combinations of hat functions Vk G /Cm- 
Furthermore, for the contour method we choose 

• I = 10 right hand sides from /C40, 

• m = 401 functional Wj = S^j S H', where S^j € M.^ is the Dirac measure 
at a;^ = -2 + j = 0, . . . . 400, 

• symmetric finite intervals J = [—■§,■§] of length L, 

• the number x for the rank test (2.15) such that the singular values of D'^' 
satisfy ai ^ . . . ^ a^t ^ Oai > a^+i for some given 6 > 0. 

There are always two eigenvalues inside the circle, the zero and the unstable eigen- 
value, see Figure 5 (a). We take the unstable eigenvalue for tests of accuracy. A 
highly accurate reference eigenfunction is computed by applying Newton's method 
to the discrete boundary eigenvalue problem (trapezoidal method on a large inter- 
val [—50,50] at small step-size Ax = 0.01). The same step-size Ax = 0.01 is used 
when solving (4.1) on equidistant grids for the contour method. 
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Figure 1. Convergence of the approximate unstable eigenvalue 
for projection bc's (+) and periodic bc's (o). 
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Figure 2. Convergence of the approximate eigenfunction for pro- 
jection boundary conditions (+) and periodic bc's (o). 



6.1. Dependence on the interval size. In our first experiment we vary the 

length of the interval for the finite boundary value problems (4.1). The other data 
is fixed. In particular, we use a large number of quadrature points (M = 100) on 
the contour and 6 = 10~* for determining >f as above to keep the influence of the 
rank test small. 

In Figure 1 we plot the distance of the approximate unstable eigenvalue A„(L) 
obtained by the contour method to the reference value A„. For the contour method 
(4.1) is solved on [—■§■,■§■] with periodic (o) and projection boundary conditions 
(+). For both boundary conditions one finds an exponential rate of convergence 
with a significantly better rate for the second one, as predicted by Theorem 4.7. 
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X 

Figure 3. First component of y„ and its approximation by the 
contour method for different interval sizes. 

The same observation is true for the convergence of the eigenfunction. This is 

shown in Figure 2, where we compare the angle between the approximate eigenfunc- 
tion yu{L) and the reference eigenfunction y^- For the approximate eigenfunction, 
we use (2.21)-(2.24) with hat functions ^^(a;) = max{0, 1 - \k - 200 - 100a;|}, 
fc = 0, 1, . . . , 400, which, together with Wk, form a biorthogonal system. 

In Figure 3 we plot the u-component of the reference eigenfunction (which is 
actually an approximation on [—50,50]) and compare it to the approximate eigen- 
functions obtained for the different interval lengths i = 5, 10, 15 with projection 
bc's. The approximate eigenfunctions are only shown on [—2, 2] which is contained 
in all domains of definition. 

6.2. Dependence on the number of quadrature points. In our second exper- 
iment we fix the interval to [—50, 50] with projection boundary conditions and vary 

the number of quadrature points. As functionals Wj wc choose the point evaluation 
at all grid points. We determine the rank as before and take 9 = 10~^° to keep 
the influence of the rank test small. All other data are the same as in the first 
experiment. The results are shown in Figure 4. 

One observes exponential convergence rate with respect to the number of quad- 
rature points, see [B] for a proof. For the eigenvalue errors there are some apparent 
resonances which have not yet been inverstigated further. It turns out that quad- 
rature errors dominate in this case and that error plots are almost identical for 
periodic boundary conditions. 

6.3. Dependence on the rank test. We now keep all data fixed but vary the 
rank test by prescribing the value of 

All other data are the same as in the previous experiment except that we choose 
45 quadrature points on the contour. In Table (c) from Figure 5 we list all singular 
values of the numerical approximation Dq in this case. In Figures 5 (a) and (b) 
we plot the approximate eigenvalues for k =2 and >f = 10, respectively. 

It turns out that the two eigenvalues inside the circle are nearly independent of 
X for X ^ 2, while the eigenvalues outside heavily depend on h. 
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Figure 4. Convergence of eigenvalue and eigenfunction with in- 
creasing number of quadrature points. 
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Figure 5. Influence of the rank test on the spectrum for two dif- 
ferent values of x (see (2.15)). o: Quadrature points, -|-: approxi- 
mate eigenvalues. 



It is shown in [B] that values outside but close to the contour still represent 
good approximations of eigenvalues. In our example, however, these eigenvalues 
are generated by the essential spectrum of the continuous problem which lies very 
close to the contour. This case is not covered by the analysis in [B] and requires 
further investigations. 
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Appendix A. Embedding results for the function spaces % and Hj 

The aim of this appendix is to show several embedding properties for the function 
spaces used in Sections 3.2 and 4. Recall the Banach spaces (3.37), (4.2), i.e. 

U = {y& L'i^, &):ye ACioc, -y' + By € L''{R, &)}. 

with norm ||?/||^ = + \\-y' + ByW^i and 

nj={yG L\J,C') : y G AC{J),-y' + By e L\J,&)}, 

with norm \\y\\l^^ = + \\ - y' + By\\%^jy 

In the following we assume B e L^(M, C*'**) and let T{x,xo), x,xo S M, denote 
the solution operator of —y' + B{x)y (that is, the propagator of the differential 
equation y' = B{x)y on R) in the sense of Caratheodory, i.e. T{-,xo) is a mild 
solution of the initial value problem Y' = B{x)Y, Y{xo) = Id in the sense of 
Caratheodory, or equivalently 

T{x, Xo) = Id+ [ B{^)T{^, xo) d^, X, xo e M. (A.l) 

Jxo 

Then T e C{R^,C'^''^) and T{-,xo) € ACioc(R)- Gronwall's inequality yields for ah 
x, Xo e R: 

\\T{x,xo)\\^exp{ [ \\B{0\\d^) ^exp{\\BU^) =:K. (A.2) 

Every y € Hj for J C R satisfies the equation —y' + By = z on J for some 
z G i^(J, C^) and hence for all a;, ccq G J C R we have 

yix) = T{x, xo)yixo) - f T{x, 0^(0 d^- (A.3) 

J Xo 

This implies for all x, xq e J C R the estimate 

\y{x) - y(.xo)| ^ K\\B\\Li{ix,xo])\y{xo 

where we use [a;,a;o] to denote the interval [a;o,a;] in case Xq < x. Indeed, (A.4) 
follows from 

\\T{x, Xo) -Id\\<. ii'||B|Ui([x,xo])> (A.5) 

and 



/" 

J xn 



T{x,Oz{^)di\\ < ^\x-xo\K\\z\\l^, (A.6) 

>Xo 

which are easily obtained from (A.l) and (A.2). 

Our first lemma shows asymptotic decay of elements in "H, that is, embedding 
oin in the space Co(R,C<') = {y G C(R,C'^) : yim.\^\^^y{x) = 0}. We supply a 
short proof and refer to [CL, Lem.3.16] for a more general case. 

Lemma A.l. Assume B G L^{R,C'^''^)- Then U C Co(R,C'^) andUj C C(J,C<^). 

Proof. Inclusions U C C(R,C'^) and Hj C C(J,C'^) follow from (A.3). Assume 

that for some y € % there is a sequence (.Tjv)jvgN C R with |xjv| — >■ oo, so that 
|2/(a;v)| ^ f > 0. Without loss of generality we may assume xat+i ^ xat + 1. Now 
let So = min{i, u^{3K\\y\\'H)~'^}. Since B G ^^(R, C'^''') we can choose < Si ^ So, 
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such that K Jj^ ||-B(a;)|| dx ^ ^ for all measurable C M with meas(AI) < 2Si. 
For |a; — a;jv| ^ inequality (A.4) then implies 

\y{x)\ > |y(a;jv)| - ^\y{xN)\ " | ^ |- 
This leads to a contradiction via the estimate 

f °° pXN + Sl 

\\y\\n^ / \y{x)\'dx^Yl / \y{x)\'dx = ^. 



Next, we show that embedding of Hj in L°°{J, C^) is uniform for all J C K. 

Lemma A. 2. There exists a constant C > such that for all intervals J = [a, b] C 
M with 6 — a > 1 and for J = 

l|y|U~(.7) ^ C\\y\\Hj, for ally€ Hj. (A.7) 

Proof. Let y G Hj, then j|t/||L=c < oo and there is a; € J with = 
by Lemma A.l. Let 6 G (0, min(l, iiT^^)) be so small that K Jj^ \B{x)\dx ^ ^ 
holds for all measurable M CR with meas(A^) < 6. Let C = 2^^{2 + 2^)/S. If 
Ibllioo ^ 8||y|l2^^ then (A.7) holds since C > y/S. If \y{x)\^ > 8\\y\\^^ then, by 
(A.4), for all x E ,J with |.x — .t| ^ (5 we have 

\yix)\ > \y{x)\ - i|y(x)| - ^rK\\y\\n, > \\y{x)\ - \\y\\n, > 0. 
Prom this we obtain 

\\y\\nj> / \y{x)fdx^ / {^\y{x)\ - WvWhj)^ dx 

J J J jn[x-s,x+s] 

> <5 f Mk _ M,^\\yu, + Ml,,) > s (ly^ - ibii^ 



and (A.7) follows. 
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